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Abstract: Parton shower Monte Carlo event generators in which the shower evolves 
from hard splittings to soft splittings generally use the leading color approximation, which 
is the leading term in an expansion in powers of where A'^c = 3 is the number of 

colors. We introduce a more general approximation, the LC+ approximation, that includes 
some of the color suppressed contributions. There is a cost: each generated event comes 
with a weight. There is a benefit: at each splitting the leading soft x collinear singularity 
and the leading collinear singularity are treated exactly with respect to color. In addition, 
an LC+ shower can start from a state of the color density matrix in which the bra state 
color and the ket state color do not match. 
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1 Introduction 

Partons carry momentum, flavor, spin, and color. All of these quantum numbers are 
represented in parton shower Monte Carlo event generators like Pythia [1], Herwig [2], 
and Sherpa [3]. Spin and color do not fit easily into the event generator format because 
quantum interference between different spin and color states is important, even in the limit 
of very soft or collinear parton splittings. In this paper, we address the issue of color. 

In previous work [4], we have shown how the evolution equations for a parton shower 
can be formulated in a way that fully includes spin and color. The resulting integrals 
can, in principle, be evaluated by Monte Carlo integration. However, simple Monte Carlo 
methods will not be practical when there are many splittings to be represented. In ref. [5], 
we found that these same evolution equations simplify if we average over spins and take the 
leading color limit, dropping all terms proportional to where A^c = 3 is the number 

of colors. Then one gets evolution equations that can be realized as a Markov process with 
positive probabilities. In ref. [6], we described a method for incorporating spin interference 
that we believe is practical. We are currently working on code to realize all of this. In 
order to keep the exposition in this paper as simple as possible, we average over spins and 
concentrate on color only. 

The problem of treating color interference is, in our experience, more difficult. One 
method is to order gluon emissions according to emission angle, which takes into account 
quite a lot of the physics of color coherence at the cost of approximating the full dependence 
of the emissions on the emission angles [7]. With the approximations, this does give positive 
splitting probabilities. The other main method is to simply drop terms that correspond to 
interference in the color space on the grounds that these terms are suppressed by factors 
This is the standard leading color (LC) approximation. 

In this paper, we describe a less restrictive approximation, which we call the LC+ 
approximation. There is a cost to using the LC+ approximation in place of the LC ap- 
proximation: one then gets shower events with weights, which can be negative. However, 
we argue that the deleterious effect of the weights on the numerical performance of the 
algorithm can be controlled - for instance by switching back to the LC approximation 
after some number of parton splittings. There are several advantages to the LC+ approx- 
imation. First, the color changes in collinear emissions and in emissions that are both soft 
and collinear are treated exactly. The only approximation is in wide angle soft emissions. 
Second, LC-|- evolution can be started from any partonic color state, including states in 
which partons in the quantum amplitude have one color configuration while partons in 
the conjugate quantum amplitude have a different color configuration. The LC approxi- 
mation cannot work with this kind of interference state. Third, the LC+ approximation 
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is quite flexible, so that it could be applied within any parton shower program based on 
color dipoles, such as Pythia or Sherpa. For this reason, we believe that the method is 
of quite general interest. 

We find that to state the method precisely, we need a fairly elaborate formalism based 
on linear operators on a certain vector space. We borrow this formalism from Refs. [4-6]. 
However, the essence of the LC+ approximation can be understood from a consideration 
of standard Feynman-like diagrams that represent color flow. For this reason, we provide 
a heuristic introduction to the approximation in sections 2, 3, and 4. We then turn to the 
more detailed, operator based, analysis in sections 5, 6, and 7. In the LC+ approach, there 



In section 8, we describe how to get back to diagonal configurations, {c'}m = {c}m at the 
end of the shower. In section 9 we describe how one can include in perturbation theory the 
effects that are neglected in the LC+ approximation. In section 10 we include the effects 
of the color phase induced by exchanging virtual soft gluons. Finally, in section 11, we 
conclude the main part of the exposition. We treat some more technical topics in three 
appendices. 

2 Statement of the problem 

In this section, we seek to describe how color evolves in a parton shower, illustrating why it 
is not simple to describe the color evolution exactly. First, let us note that, in a perturbative 
treatment, a parton carries momentum, flavor, spin and color. We can imagine that the 
momenta of partons at the end of the shower are quite precisely measured. Then the 
momenta of intermediate partons in a shower are also well determined. Thus we do not 
need to consider interference between states in which intermediate partons have different 
momenta p and p' . Since the flavor of a mother parton in a splitting is determined by 
the flavors of its daughters, we also do not need to consider interference between states 
in which intermediate partons have different flavors / and /'. However, this argument 
does not hold for spin and color. For instance, in a splitting 1—7-2 + 3 for gluons with 
colors ai, 02 and 03, the matrix element is proportional to /ai ,02,03 and for fixed 02 and 03, 
/ai ,02,03 can be nonzero for more than one value of ai. Thus we need to allow for quantum 
interference between different color states of the same parton. We also need to allow for 
quantum interference between different spin states. 

2.1 Evolution without color or spin 

Let us consider what happens in a parton shower that evolves from hard splittings to 
soft splittings.^ To get started, we ignore both spin and color. We define a "shower 
time" variable t such that an initial hard parton scattering happens at t = and then at 
each interval dt a parton has some probability of splitting to become two partons. Harder 
splittings happen at smaller t values, successively softer splittings happen at larger t values. 
For instance, exp(— t) can be proportional to the virtuality in a splitting or to the transverse 

^We have in mind something like Pythia or Sherpa. The program Herwig, based on ordering in angles, 
is rather different. 



are different color states for the amplitude 
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momentum of the daughters relative to the mother parton direction. As implemented in 
a computer program, the partonic system always has a definite state. Ignoring spin and 
color, the state of the system is described by the momenta and flavors. 



{PJ} 



m — 



{Pa, fa.,Ph, fh,Pl, fl, ■ ■ ■ ,Pm, fm} ■ 



(2.1) 



Here Pa is the momentum of the incoming parton from hadron A, equal to a fraction rjg^ of 
the hadron momentum, pb is the momentum of the incoming parton from hadron B, and 
the Pi are the momenta of m outgoing partons. The flavors of the partons are denoted by 
discrete flavor variables /. In each shower time interval dt, there is a certain probability 
that the system will switch to a new state. What actually happens in a given simulated 
event is determined by generated pseudo-random numbers. This means that in an ensemble 
of simulated events, there is an probability pi{p, f}m,t) that the partonic system is in a 
certain state at shower time t. This probability distribution then evolves with t. Let us 
denote by \p{t)) the function p at time t considered as a vector in the space of functions of 
m and {p, f}m- Then the evolution equation for the probabilities has the form of a linear 
equation 



Here Tiiit) is a linear operator on the space of probability distributions. The operator 
'Hi{t) corresponds to the splitting probabilities chosen for shower evolution. (There are 
many possible choices.) Then V{t) is another linear operator that is constructed from T-Li{t) 
so as to conserve probability in the shower evolution. The Sudakov factor that represents 
the probability that there was no splitting between times ti and t2 is exp(— j^^ dt V{t)). 
We will return to this in later sections. For the moment, all that we need to know is that 
the evolution of the probability function p is determined by parton splitting probabilities. 

2.2 Ignoring spin 

Now we need to consider spin and color. We have addressed the spin problem in ref. [6]. 
Combining the spin treatment of ref. [6] with the discussion of this paper is not difficult, 
but adds a layer of complexity. In order to keep this paper as simple as possible, we ignore 
spin by supposing that we sum over the spins of the daughter partons in a splitting and 
average over the spins of the mother parton. Thus in this paper we address the color 
problem but not the spin problem. 

2.3 The color density operator 

The natural language for describing an evolving probability distribution is statistical me- 
chanics. Indeed, eq. (2.2) is a standard sort of equation for evolution in statistical me- 
chanics. In this paper, we want to include quantum color, so we need quantum statistical 
mechanics. Thus we need a density operator that depends on t and is an operator on the 
space of color states of a possibly large number of partons. The density operator then has 
the form 






(2.3) 



m 
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Figure 1. Illustration of an open string color basis state (left) and a closed sting color basis state 
(right). These diagrams represent the color matrices without the normalization factors N(S)^^^'^ . 

Here |{c}m) and |{c'}m) are basis vectors for the quantum color space. The color configu- 
ration {c}m of the quantum amplitude is, in general, different from the color configuration 
{c'}m of the conjugate quantum amplitude. Thus the function p{{p, f,c' ,c}m,t) depends 
on two sets of color indices. The density operator p({p, f}m,'t) can be regarded as a vector 
in the space of functions of {p, f}m with values in the space of operators on the quantum 
color space. It obeys an evolution equation of the form (2.2). We simply have to reinterpret 
what \p{t)) means. Of course, with color involved, the detailed structure of eq. (2.2) is 
now more complicated. To understand the color structure, we need to first say something 
about the color basis states. 

Our analysis in this paper makes use of the standard color basis defined by two sorts 
of vectors [4, 8]. First, there are open string vectors 

^(5)^"> = NiSy^/^i^H"^ ■ ■ ■ i""-^]ai,a„ • (2.4) 

Here oi is a quark color index, a„ is an antiquark color index, and the other are gluon 
color indices. The are standard SU(3) color matrices in the fundamental representation. 
Also, N{S) = NcCp~'^ is a normalization factor such that (^|^) = 1- Second, there are 
closed string vectors 

^'(5){4 = 7v(5)-i/2xr[t"it'^2 ...t'^"] . (2.5) 

Here all of the are gluon color indices. Also, N{S) = Cp is a normalization factor such 
that (^1^) = 1 — [~^/{^c ~ This normalization factor is close to 1, with a small 

correction. The most general color basis vector, which we denote by |{c}m), is a product 
of these two kinds of units. The open string and closed string basis vectors are illustrated 
in figure 1. 

The color basis states are normalized to ({c}m|{c}m) = 1 or to ({c}m|{c}m) ~ 1 with 
very small corrections. They are not, however, generally orthogonal. However, when {c}m 
and {c'}m are different, one finds that ({c'}m|{c}m) = 0{l/N^). That is, the basis vectors 
are orthogonal in the Nc ^ oo limit. 

Now let's look at an example. In figure 2, we depict |{c}m)({c'}m| for a color structure 
that arises in the hard scattering process q{a)+q{h) — )• g{l)+g{2). Both |{c}m) and ({c'}m| 
show the color state that corresponds to t-channel q exchange. Thus we have a diagonal 
contribution to the color density operator, with {c}m = {c'}m- 
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|{c}m) {{c'}m\ 

Figure 2. Illustration of the color density operator. We depict |{c}to)({c'}„i| with a quark, an 
antiquark, and two gluons. Here {c}m — {c'}m- 




Figure 3. Illustration of the color density operator as in figure 2. In the case illustrated, we have 
{c}m 7^ {c'}m because we have switched the positions of gluons 1 and 2 in {c'}m- 




a 





Figure 4. Identity for color dependence of gluon splitting. 



One can also have li-channel q exchange, which amounts to exchanging gluons 1 and 
2. In figure 3, we still have the t-channel diagram for {c}m, but now we illustrate the 
contribution from the u-channel diagram for {c'}m. This gives an ofF-diagonal contribution 
to the color density operator, with {c}m / {c'}m- Each of the two contributions to the 
color density operator shown in the two figures can serve as the starting point for a parton 
shower. Since their color structures are different, they should generate different showering. 
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Figure 6. Identity for color dependence of splitting of gluon 1 in the ket state with the participation 
of helper parton b in the bra state. 

2.4 Color structure of the parton splitting operator 

Having seen the meaning of the color density operator, we can now consider what happens 
in shower evolution when a gluon splits. In figure 4, we show the color structure when 
mother gluon 1 splits into daughter gluon 1 and daughter gluon 3. Using the identity 
ifaia^c tc = ^ai^as " ^as^ai > wc find that there is a term in which the new gluon 3 is inserted 
below gluon 1 on the quark line minus a term in which the new gluon is inserted above 
gluon 1 on the quark line.^ 

Applying this identity, we see in figure 5 that letting mother gluon 1 split in both the 
ket state and the bra state in figure 2 gives contributions with four new |{c}m+i)({c'}m+i| 
terms. Even when we use only the simplest kind of splitting and even though we start 
with a color diagonal density operator, we generate off diagonal terms in the color density 
operator. 

There can also be quantum interference between emissions of a soft gluon from two 

^In our diagrams, the three gluon vertex is ifaiasc with (aiOac) ordered clockwise around the vertex. 
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Figure 7. Identity for color dependence of splitting of gluon 1 in the ket state with the participation 
of helper parton 2 in the bra state. 

different partons. We illustrate this in figure 6. Mother parton 1 emits gluon 3 in the ket 
state, while the antiquark with label b emits gluon 3 in the bra state. In a physical gauge, 
the matrix elements for this are large as long as gluon 3 is soft. One can consider that a 
color dipole consisting of partons 1 and b emits the gluon. In a dipole shower, one wants 
to account for emission from each dipole (insofar as complications from color allow). In 
a dipole antenna shower like ViNCiA [9], one keeps the dipoles together as a unit. In this 
paper, we follow the method of a partitioned dipole shower like Pythia, which partitions 
the radiation pattern into two terms. One term is singular when the soft gluon is collinear 
with parton 1, the emitting parton, but is not singular when the soft gluon is collinear with 
parton b, the helper parton. The other term is singular when the soft gluon is collinear 
with parton b but is not singular when the soft gluon is collinear with parton 1. For our 
present discussion, let us consider parton 1 to be the emitting parton. On the right hand 
side of the figure, we have applied the color identity of figure 4 to expand this contribution 
in our standard color basis. We see that we get two off diagonal contributions. 

Parton 3 can also be emitted from parton 1 with parton 2 as the helper parton. This 
case is illustrated in figure 7. When we expand in five-parton color basis states, there are 
four contributions. One of them is color diagonal. 

Let us summarize. Using a few simple examples, we have explored the color content 
of the splitting operator Hi{t) in eq. (2.2). There are two sorts of graph, direct graphs 
like figure 5 and interference graphs like figures 6 and 7. The color structure of both kinds 
of graphs is non-trivial. Even if we start from a color diagonal contribution to the color 
density operator, each splitting creates off-diagonal contributions. With each new splitting, 
the color state becomes more complicated. 
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2.5 Color structure of the virtual splitting operator 

What about the operator V{t) in eq. (2.2)? This operator maps a state of the color 
density operator with m final state partons into another state with m final state partons. 
Physically, it represents virtual Feynman graphs with hardness scale characteristic of the 
current shower time t. We determine what V{t) should be by demanding that shower 
evolution conserve probability. This means that 

{l\[niit)-V{t)]\pit))=0 , (2.6) 

where multiplying by (l| is a convenient way of writing the instruction to take the trace 
of the color statistical operator and integrate over all of the parton variables. 
Note that taking the trace of the color statistical operator means using 

Ti[p{{p,f},n,t)]= Yl Pi{pJ,c,c'U,t){{c'y^\{cU) . (2.7) 

{c}m,{c'}m 

We get the inner product ({c'}m|{c}m). It seems unfortunate that ({c'}m|{c}m) does not 
vanish when {c'}m 7^ {c}m- However, this inner product is suppressed by powers of l/N^ 
when {c'}m / {c}m- 

It is possible to satisfy eq. (2.7) and, at the same time, match the color structure of 
virtual Feynman graphs by letting V have the form 

V = (/i + i(/>) ® 1 + 1 (/it - i^) . (2.8) 

This notation denotes that {h + 10) is an operator on the ket part of the color density 
operator and {h^ — icj)) is an operator on the bra part. The operator (j) is hermitian and 
represents a color phase. We will consider in section 10. Until then, we simply set = 0. 
The operator h is hermitian in the full color theory, but the LC+ approximation for it is 
not. For that reason, we distinguish the roles of h and in our formulas. With = 0, we 
have 

V = /i (g) 1 + 1 /it . (2.9) 

We can state this a little more precisely. If p{{p, f}m,t) is defined by eq. (2.3) and \p') = 
V\p), then 

P{{pj}m,t)= ^ p{{p,f,c,c}m,t) 

{c}m,W}m (2.10) 

X [/l(fe/}m)|{c}™>({c'}™| + |{c},„>({cV|/lt({p,/}„)' . 

Each term in Tij determines a contribution to h. To see how this works, consider the 
splitting in figure 5. We define h so that the corresponding contribution to has the 
color factor shown in figure 8. To verify probability conservation, eq. (2.6), just take the 
color trace of figure 8 and compare it to the color trace of the left hand side of figure 5. In 
each case we get the number represented by the color diagram in figure 9. 

The color structure of V is more complicated in the case of interference diagrams. 
Consider the splitting in figure 7, in which parton 1 in the ket state splits and we take the 
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Figure 8. Color structure of contribution to V{t) from the splitting in figure 5. 




Figure 9. The color trace of figure 8 of the left hand side of figure 5. 




Figure 10. Color structure of contribution to V{t) from the splitting in figure 7. 



interference with the emission of the same gluon from helper parton 2 in the bra state. We 
can consider this diagram to correspond to a contribution to V in which h acts on the bra 
state. We define this contribution to h so that jp') has the color factor shown in figure 10. 
To verify probability conservation, just take the color trace of figure 10 and compare it to 
the color trace of the left hand side of figure 7. Then the contribution to V in which h acts 
on the ket state corresponds to the complex conjugate of the splitting in figure 7. 

Note that the contribution to V{t) corresponding to a direct splitting term in is 
trivial. Evidently from figure 8, this contribution to V(t) simply multiplies the color state 
by Ca- In contrast, the contribution to V{t) corresponding to an interference term in Hj 
is not trivial. When one expands the bra state in figure 10 in the color basis states, we 
get several contributions.^ Thus, in general, V{t) is not diagonal in color. This makes an 
exact accounting for color difficult because V{t) enters the generation of parton showers in 
the form of the Sudakov factor, the time ordered exponential Texp — dr V(r) 



^For the case depicted, there are two contributions. In other cases, there are more. 
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3 The leading color approximation 

We have now seen something of the structure in color space of a leading order parton 
shower (of the partitioned dipole variety). Of course, this is not the structure of any 
existing computer program that models a parton shower with quantum color. There is no 
such program. The evolution equations make good sense, but implementing more than a 
couple of splitting steps as a computation is not feasible with current methods. There is, 
however, a simple approximation that one can apply to get a practical implementation. 
This is the leading color (LC) approximation. 

The general idea of the leading color approximation is to neglect contributions to the 
probability to get a given final partonic state that are suppressed by powers of l/N^- 
To proceed, one replaces gluons in the 8 representation of SU{3) by gluons in the 3x3 
representation, using 

tijt°;,j, - 6ij' 6i'j . (3.1) 

One also notes that all contributions to the color density operator with {c}m 7^ {c'}m 
can simply be dropped. That is because ({c'}m|{c}m) is suppressed by a power of l/N^ 
compared to ({c}m|{c}m) or {{c'}m\{c'}m)- During shower evolution, the ket states and 
the bra states evolve separately, but at the end of the parton shower, when there are 
final state partons, one needs to use ({c'}iv|{c}iv) to compute a probability. With a little 
analysis (see section 5.11), one sees that once one has {c}m W}m, one can never get 
back to a leading power by further parton splittings. The consequence of this is that 

as splitting proceeds, one simply drops {c}m 7^ {c'}m contributions. Thus on the right 
hand side of figure 5, one keeps the first two terms and drops the second two terms. On 
the right hand side of figure 7, one keeps the first term and drops the other three. On the 
right hand side of figure 6, one drops both terms. That is, interference between emission 
of a gluon from parton / and from a helper parton k that is not directly color connected 
to I is neglected. 

With the leading color approximation, the no-splitting operator V is simple. For 
instance, on the right hand side of figure 5 we keep the first two terms. In each of these 
terms, when we take the color trace we see that adding the emitted gluon 3 creates one 
more color loop and thus one more factor of Nc- Thus we can take the color operator in 
V to simply multiply the color state we started with by Ca = Nc (including two graphs, 
each with a factor 1/2 from eq. (3.1)). 

For a splitting q — t- q + g, the color factor in V would be Nc/2 when we calculate this 
way, but one normally uses Cp instead by simply multiplying the splitting probability by 
2Cf/Ca. 

Strictly speaking, one would drop all g ^ q + q splittings in the leading color approxi- 
mation, but one normally keeps the part of this splitting suppressed by i/Nc, omitting the 
parts suppressed by l/N^. 

The leading color approximation gives a simple shower algorithm. It is also intuitively 
appealing. The ket states and the bra states in the color density operator always have the 
same {c}m- In {cjm, we can think of each quark or antiquark as being connected by a color 
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string to a neighboring gluon in a color basis state in figure 1. Each gluon is connected 
to two other partons by a color string. Then emitting a new gluon means connecting it 
to two of the previous strings. Interference between emitting a new gluon from parton / 
and helper parton k can only occur if the two partons were color connected; then the new 
gluon is connected to the string that previously joined partons / and k. 

One can note two problems with the leading color approximation. First, it neglects 
terms suppressed by powers of Second, it cannot start with color density operator 

contributions |{c}m)({c'}m| with {c}m / {c'}m- 

4 Introduction to the LC+ approximation 

We propose in this section an improved "LC+" approximation that goes beyond the lead- 
ing color approximation by including some of the contributions to cross sections that are 
suppressed by powers of 

To go beyond the leading color approximation, one has to give up something.^ We 
choose to give up having an algorithm that can be implemented without having weights 
for events. In particular, the terms in the splitting probabilities in the LC+ approximation 
have both plus and minus signs. One cannot generate events with negative probabilities, 
so in order to include these terms one will have to generate certain splittings with positive 
probabilities and negative weights. The weights are then carried with the event. Now, 
having weighted events is not intrinsically a problem for convergence of Monte Carlo inte- 
grations to a physical answer. However, there would be a problem that would not be easy 
to avoid if a physical cross section to be calculated had the form a"o[1.00 — 9.01/N^ + •••]. 
Then calculating the result by Monte Carlo integration would not give an accurate answer 
for A'c = 3. We believe that this does not happen. If it does happen for some observable 
cross section, then the LC+ approximation may not work well. Of course, in this case, the 
standard LC approximation is giving us the wrong answer and the LC+ result would at 
least give an indication of problems. 

We will state what the LC+ approximation is quite precisely using an operator notation 
in section 6. However the main idea can be grasped quite easily using the examples that 
we have already explored, so we do that in this section. 

4.1 The parton splitting operator 

First, the direct graph for the splitting of parton I by emitting a gluon labeled m + 1 
contains (in a physical gauge) the collinear singularity that occurs when the parton mo- 
menta after the splitting obey Pm+i pi and also the double singularity that occurs when 
Pm+i — ^ with Pm+i oc pi- We want the coefficients of the logarithms that arise from these 
singularities to be exactly right with respect to color. Thus we keep the color structure of 
these splittings exactly. This means that in the right hand side of figure 5 we keep all four 
terms. We can start with {c'}m 7^ {c}m as in figure 11. Again, there are four terms and 
we keep them all. 

*More precisely, the authors do not know how to go beyond the leading color approximation while giving 
up nothing. 
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Figure 11. Identity for color dependence of splitting of gluon 1 in both the bra state and the ket 
state in the case {c'}m 7^ {c}m- In the LC+ approximation, one keeps all four terms. 




Figure 12. LC+ approximation for splitting of gluon 1 in the ket state with the participation of 
helper parton 2 in the bra state in a case with {c'}m 7^ {c}m- 

Now consider a splitting in which gluon emission from parton I in the ket state interferes 
with emission from helper parton k in the bra state. An example was shown in figure 6. 
Here the helper parton (A; = b) is not color connected to the emitting parton (/ = 1) in 
the bra state. In this case, the LC+ approximation is to drop the contribution entirely. A 
second example was shown in figure 7. Here the helper parton {k = 2) is color connected 
to the emitting parton (/ = 1) in the bra state. In this case, the LC+ approximation is 
to keep the contributions in which the emitted gluon attaches between / and k in the bra 
state. That is, we keep the first and third terms on the right hand side of figure 7. A third 
example is shown in figure 12. In this case, {c'}m 7^ {c}m- The helper parton (k = 2) is 
color connected to the emitting parton (/ = 1) in the bra state, so the LC+ approximation 
is to keep the contributions in which the emitted gluon attaches between I and k in the 
bra state. The two contributions that are retained in the LC+ approximation are shown 
on the right hand side of figure 12. 
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Figure 13. Color structure of contribution to V{t) from the splitting in figure 11. 




Figure 14. LC+ approximation for splitting of gluon 1 in the ket state with the participation of 
helper parton 2 in the bra state in a case with {c'}m 7^ {c}m- This is the sum of the two terms on 
the right hand side of figure 12. 




Figure 15. Color structure of contribution to V{t) from the splitting in figure 14. 
4.2 The virtual splitting operator 

To find the structure of the virtual splitting operator V{t) in the LC+ approximation 
we use the definition (2.9). For a direct splitting diagram, we consider the real emission 
diagram to correspond to the sum of virtual diagrams in which h acts on the ket state and 
in which /i^ acts on the bra state. For an interference diagram with the helper parton in 
the bra state, we consider the real emission diagram to correspond to a virtual diagram 
in which h acts on the bra state. An interference diagram with the helper parton in the 
ket state corresponds to a virtual diagram in which h acts on the ket state. Then we 
impose probability conservation, eq. (2.6), to relate V(t) to the definitions of the LC+ 
approximation for 'Hi{t) outlined above. 

For a direct splitting diagram as in figure 11, this gives the contribution to V{t) illus- 
trated in figure 13. To see this, we simply note that the color trace of the left hand side of 
figure 11 is the same as the color trace of figure 13. 
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Figure 16. LC+ approximation for splitting of parton "a" in the ket state with the participation 
of helper parton 2 in the bra state in a case with {c'}m ^ {c}m- 

For the interference diagram in figure 12, we first note that the two terms that are 
retained in the LC+ approximation sum to the contribution shown in figure 14. This gives 
the contribution to V(t) illustrated in figure 15 since the color trace of the right hand side 
of figure 14 is the same as the color trace of figure 15. 

We now note something remarkable. The color basis states are eigenstates of the 
operator h that defines V{t). For a direct splitting of a gluon, as in figure 13, the eigenvalue 
is Ca/2. Similarly, for a direct splitting of a quark or antiquark, the eigenvalue is Cf/2. 
For an interference diagram in which a gluon splits to two gluons with interference from 
gluon emission some other parton, as in figure 15, a simple calculation shows that the 
eigenvalue is Ca/4. Similarly, for an interference diagram in which a quark or antiquark 
splits, emitting a gluon with interference from gluon emission some other parton, the LC+ 
approximation is shown in figure 16. Then the corresponding contribution to V{t) has 
eigenvalue Cf/2. 

5 Operator based analysis with full color 

In the preceding sections, we have sketched the main idea of the LC+ approximation. Now 
we need to make the idea precise. To do that, we first set up a formalism to describe 
the evolution of generic parton shower that includes quantum color exactly. To keep the 
presentation simple, we average over spins. More precisely, we sum over the spins of partons 
after each splitting and average over the spins of the partons that enter each splitting. Then 
spin is not visible in the evolution equations at all. The formalism is taken from Refs. [4- 
6]. However, here we keep the shower quite generic by not specifying exactly the splitting 
functions, the definition of the shower time t that is used to order successive splittings, or 
the momentum mapping for each splitting that allows all partons to be on-shell while at 
the same time exactly conserving momentum. 

5.1 Parton labels 

At each stage of the shower, there are two initial state partons with labels "a" and "b" 
together with m final state partons with labels 1, . . . , m. The parton momenta and flavors 
are then specified by a list {p,f}m as in eq. (2.1).^ At each step of the parton shower, 

^For the flavors for the initial state partons, it is useful to let / denote the flavor leaving the hard 
interaction, which is the opposite of the flavor of the physical parton entering the hard interaction. 
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any one of the partons can split. This includes an initial state parton, which splits in 
"backwards evolution" to another initial state parton plus a final state parton. In either 
case, let I be the label of the parton that splits. At a splitting, the parton / remains and 
one more parton, with label m + 1 is created. After the splitting, the momenta and flavors 
are {p,f}m+i- Whenever a final state gluon is created, we assign the label m + 1 to the 
gluon. In the case of a final state g ^ g + g splitting, one daughter gluon has label / and 
the other has label m + 1. We use the interchange symmetry of the process to rearrange 
the splitting function so that there is a singularity when gluon m + 1 becomes soft but no 
singularity when gluon I becomes soft. 

In the case of gluon emission, there are interference graphs. A gluon m + 1 emitted 
from parton / in the ket state can be emitted by partner parton k in the bra state. Similarly 
a gluon m + 1 emitted from parton I in the bra state can be emitted by partner parton k 
in the ket state. The interference diagrams are important when gluon m + 1 is soft. 

5.2 The evolution equation 

The development of a parton shower can be described by giving an evolution equation for 
the operator p{{p, f}rmt) defined by eq. (2.3) as the density operator on the color space 
of the evolving partonic system. (It would be an operator on the spin space also except 
that we average over spins in this paper.) The density operator is a function of the parton 
momenta and flavors. This function can be regarded as a vector \p{t)) in the space of 
functions of {p, f}m with values in the space of operators on the quantum color space. We 
take the vector |/o(t)) to obey a linear evolution equation of the form (2.2). 

In eq. (2.2), 'Hi{t) is the splitting operator. Specifying shower evolution means speci- 
fying what 'Hi{t) is. To do that, we use basis states /, c', c}m) that specify a definite 
momentum, flavor, and color configuration for the partons. The normalizations are such 
that p({p, /,c',c}m,t) in eq. (2.3) is 

P({P, /,c',c}m,t) = ({p,/,c',c}m|/j(t)) . (5.1) 

Making use of these basis vectors, we can define T-Li{t) by specifying the matrix elements 
of Tii{t) between a partonic state after a splitting and the state before the splitting, 
({p, /, c', c}m+i|'H/(f)|{p, /, c', c}m) . We begin this task in the next few subsections by 
discussing some of the ingredients in this matrix element. 

5.3 Splitting functions 

To describe emission of a soft gluon m + 1 from parton I with interference from emission 
from helper parton k, we use the dipole splitting function 



— Qipoie . -{Pm+1 -piPk- Pm+1 ' Pk Pi? ,^ „s 

^ik = 47ras ^ ^rT2 • (5.2) 

[Pm+l ■ Pk Pm+1 ■ Pi) 



- dipole 



In this expression, partons I and k can have nonzero masses. The expression for wik may 
be more familiar in the massless limit, pf ^ and p| — t- 0, where it becomes 

— dipole , . "^Pk 'Pi r,\ 

W]i, 47rQs ; — . (5.3) 

Pm+1 ■ Pk Pm+1 ■ Pi 
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Eq. (5.2) includes all four diagrams for emission from either parton / or parton k in both 
the bra and ket states, calculated in the limit Pm+i = ^Pm+i with A — >• 0. If we calculate 
the four contributions separately using the eikonal approximation in a physical gauge, then 

-dipole ^ -eikonal ^ ^^ikonal _ 2^ikonal (5 4) 

The first two terms represent the direct graphs while the last term represents the two 
interference graphs.^ 

The splitting function ^tj^ikonal (5.4) is singular in the limit in which parton I is 

massless and parton m + 1 is collinear with parton /: "Pm+i [(1 — ^)/APl- this limit, 

eikonal ^ 8vr«s 2z 

2p^+i -pil-z 

There is a factor of the virtuality of the splitting, 2pm^i -pi in the denominator and there is 
a function of the momentum fraction z that is singular in the soft gluon limit z — )■ 1. The 
behavior of this function away from the collinear limit depends on the conventions used to 
define the shower. For instance, it depends on the definition of z. 

The splitting of parton / into a parton of the same flavor plus a gluon both in the ket 
state and in the bra state is described by a function wu. Suppose that parton / is a quark. 
Then in the limit in which parton / is massless and parton m + 1 is collinear with parton 
I we have 

_ StTOs l + z^ 

wii^— . 5.6 

2pm+i -pi I - z 

Away from this limit, the exact result depends on the conventions used to define the 
shower.^ Note that matches wu in the limit of a collinear splitting that is also soft, 

(l-z)^O. 

For other sorts of splittings, there are splitting functions wu that are proportional to 
the DGLAP splitting kernels in the limit of massless, collinear splittings. 

5.4 Shower time 

Parton showers are based on evolution of the system as a shower time variable t increases. 
The idea is to start at the hard interaction and move to a first splitting that is less hard, 
then move to softer splittings as t increases. For initial state splittings, this means moving 
backward in physical time. 

As a measure of softness, one can use the virtuality of the splitting, the virtuality 
divided by the energy of the mother parton, or the square of the part of the momentum 
of one of the daughters that is transverse to the direction of the mother. The shower 



®The function wf^°''''\ which equals «;fe*°"''', is denoted by Wik/{2Aik) and given in eq. (4.4) of ref. [5]. 
The function wf;'^"""' is denoted by Wa ""'^ and given in eq. (2.10) of ref. [5]. These functions inchide only 
the leading contribution in the limit A — > 0. One could, in principle, add contributions proportional to 
A^^ and higher powers of A. 

^For a g ^ g + g splitting, the DGLAP splitting kernel appears in wu, as in eq. (5.6), once one 
symmetrizes in z 1 — z. Before symmetrization, the result depends on the conventions used to define the 
shower. 
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program Herwig orders splittings according to the angle between the daughters times the 
energy of the mother, with wide angle splittings first. This ordering has its advantages, but 
the generic shower scheme outlined here would need some modification to fit with angular 
ordering. 

5.5 Momentum mapping 

At each splitting, one starts with m final state partons with momenta {p}m and ends with 
m + 1 final state partons with momenta {p}rn+i- One might like to define pj = pj for j / I 
and j m + 1 and to take pi + Pm+i = Pi (or pi — Pm+i = Pi for an initial state splitting). 
However, this would not allow all three of pi, Pm+i, and pi to be on-shell. Accordingly, 
one needs to take some momentum from the partons other than I in order to supply the 
needed momentum in the splitting and keep exact momentum conservation. Thus we need a 
momentum mapping in which the {p}m plus three splitting variables determine the {j5}m+i- 
The three splitting variables can be, for instance, the virtuality {pi ±pm+i)^ — fn^ifl), a 
momentum fraction z and an azimuthal angle of the daughters about the direction of the 
mother. Equivalently, the {p}m+i determine three splitting variables and the {p}m- The 
needed momentum mapping can be specified by giving a function ({p, /}m+i|'P/|{P) f}m) 
that consists of an integration over the three chosen splitting variables and a product 
of delta functions that determine the {p}m+i from the splitting variables and the {p}m- 
There are many ways to do this, but for our purposes, all we need to know is that defining 
a shower evolution entails the specification of ({p, f}m+i\'Pi\{p, f}m)- 

For our generic shower, we take the mapping operator Vi to depend on the label I of 
the parton that splits. One can also let the mapping operator depend on the index k of 
the helping parton involved in interference diagrams. This is a common choice because it 
allows all of the momentum transfer to come from changing the momentum of parton k. 
That is, one can take pi + pm+i — Pi = Pk — Pk- However, this scheme is a bit awkward for 
cases in which there is no helper parton k, such as 5 — >■ g + g. Accordingly, we keep the 
generic shower equations simple by letting Vi depend only on the index I. 
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5.6 Splitting operator 

Now we are ready to specify the matrix elements of T-Li{i): 

{{V, f, c', c}m+i \'Hi{t)\{p, f, c', c}„) 

= Y,S{t- T{{p, fU+i)) {m + 1) {{p, fU+i \ri\{p, flm) 

l,k 

nc{a)nc{b) rj^rj-t, fa/AiVa., lJ^F)fb/B^'^'' I^f) 



nc(a)nc(6) 7}a57b /a/A(??a, fJ'F)fb/B{Vh, fJ'p) 

1 



X - 
2 



= ^(/m+1 / 9) Wll{{p, f}m+l) 

(5-7) 



{{C,c}m+l\t]{fl fl + /m+l) ^tkifk fk + fm+l)\{c,c}m) 
+ {{c,c}m+l\tl{fk fk + frn+l) ^ U{fl fl + /m+l) | {c , c}m) 

In the first line on the right hand side of eq. (5.7), we have a sum over indices / and k 
of partons. The parton with label / is the one that splits. There is another label k so that 
we can include graphs that represent quantum interference between emission of a gluon 
from parton I and from another parton k. We call parton k the helper parton. For the 
quantum interference terms, we have k ^ I. There are also graphs that do not represent 
quantum interference. For these, k = I. 

The first line on the right hand side of eq. (5.7) also includes a delta function that 
specifies the definition of the shower time t, then the function that defines the momentum 
mapping. 

In the second line, we have a ratio of parton distribution functions. For a final state 
splitting, this ratio is 1. For an initial state splitting, this ratio replaces the parton distri- 
bution functions at the previous momentum fraction r].^ or ?7b by the parton distribution 
function at the new momentum fraction after the splitting. 

In the following lines of eq. (5.7) there are three terms, corresponding to three types 
of splittings. 

First, there are splittings in which parton m + 1 is not a gluon. For example, a final 
state g ^ q + q splitting is in this class. For these cases, there is a splitting function 
^«/({P) /}m+i) that is singular for a collinear massless splitting, but does not have a soft 
parton singularity. 

In the second term in eq. (5.7), there are splittings in which parton I emits a gluon. The 
splitting function here, wu — w'if^°^^'^, is singular for a collinear splitting, but the singularity 
when the emitted gluon is soft has been removed by the subtraction. 

The most interesting term in eq. (5.7) is the third, in which a soft gluon m+1 is emitted 
from a dipole consisting of partons I and k. The splitting function w"^^"^^ is calculated using 
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the eikonal approximation according to eq. (5.2). It is singular when the gluon is soft and 
also contains singularities when the direction of the gluon momentum is collinear to either 
the direction of pi (if pf = 0) or the direction of pk (if pj, = 0). It is symmetric under 
interchange of pi and . 

We have manipulated the third term in order to separate the roles of partons k and /. 
Let be a function of the momenta {p}rn+i and A'^.^ be the same function with the roles 
of pk and pi reversed. Furthermore, let + ^'^^ = 1. Then, by interchanging the names 
of the dummy indices I and fc, we have 



ik^ik 



dipole 



k^l 



This gives the factor A'^f, that multiplies wfl^^°^^ in eq. (5.7). If we were to take A'^^ = 1/2, 
this manipulation would do nothing, but instead we choose A'^f^ so that = 1 when 
Pm+i is parallel to pi and = when Pm+i is parallel to pk- Then A'^jufl^^"^'^ is singular 
when the gluon is soft or collinear with parton / but not when it is collinear with parton 
k. Our preferred choice for A'^^ is given in eq. (7.12) of ref. [6]. One can also use the 
Catani-Seymour dipole partitioning function [10] for this purpose. 

Finally in eq. (5.7) there is a color factor, which is the main focus of this paper. There 
are two terms in the color factor, which are related by interchanging the indices I and k. If 
k = I, the two terms are identical. In the first term, there is an operator tj{fi — )• fi + fm+i) 
that acts on color ket states |{c}m)- This operator attaches the proper color matrix for the 
splitting to the color index for parton I. Similarly, we apply the operator tk{fk fk+fm+i) 
to the bra color state ({c'}m|. The resulting color bras and kets can be expanded in color 
basis states: 

tjifl ^ fl + fni+l)\{c}m){{c}m\tk{fk ^ fk + fm+l) 

= Yl Ci{c',cU+i,{c,cU)\{cU+i){{c'U+i\ . (5-9) 

{c',c}m + l 

Then the matrix element that we want is the corresponding expansion coefficient: 

{{C,c}m+l\t]{fl ^ fl + fm+l) ®tk{fk ^fk + fm+l)\{c,c}m) 



(5-10) 

C{{c,c}m+l,{c',c}m) ■ 



5.7 The virtual splitting operator 



The virtual splitting operator V(t) leaves the number of partons unchanged and does not 
change their momenta or flavors. It does, however, multiply by a matrix in color space. 
This matrix has two terms. First, there can be a correction to the ket color state with no 
change to the bra colors. Then, there can be a correction to the bra color state with no 
change to the ket colors. Thus the color structure of V{t) is the color structure of one loop 
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virtual corrections. We write 



V(.t)\{P,f,C,c}m) = Hl{{p, f}m; {c}m, {c}m;t)\{p, f, C , C}m) 

+ 2_/ ^rUP^ f}m;{c}m,{c}m;t)\{p, f,C ,c}m) ■ 

{&'}m 

The color matrices Hi and Hfi are the matrices that represent operators h + i4> and 
[h + = — icj) that act on the ket color vectors and the bra color vectors, respectively: 

[H{pJ}m,t)+i(l){{p,f}m,t)]\{c}m)= ^ Hl{{p, Dm] {c}ni, {c}m;t) , 

{c}m 

{{c}m\[hH{p,fU,t)-mP,fU,t)]= J2 HR{{pJ}m;{c'U,{c}m;t) {{c'U\ . 

(5.12) 

Here (/> = (f)^ . In the simplest formulation, (f> = 0. We consider (/) ^ in section 10. In 
the full color treatment, we will have = h. However in the LC+ approximation we will 
define h^^+ with [/i^^+]t ^ h^c+ . 

It is useful to solve the shower evolution equation (2.2) in the form 



\p{t)) =U{t,to)\p{to)) , (5.13) 



where 



U{t,to)=M{t,to)+ [ dTU{t,T)ni{T)M{T,to) . (5.14) 
Jto 

Here M{t2,ti) is the no-splitting operator, 



N{t2M) = Texp 



rt2 

/ dT V(t) 
Jti 



(5.15) 



that provides the Sudakov factor that is usually interpreted as the probability to not have a 
splitting between time ti and time t2- In J\f{t2, ti), the virtual splitting operators V(r) are 
time ordered. Iterating eq. (5.14) gives |p(t)) as a power series in Tij, with Tij evaluated 
at times U and factors of A/'(ii+i,ij) in between. We interpret eq. (5.14) as saying that 
possibly the system gets from to to t without splitting. If not, it goes from to to r without 
splitting, then splits according to Hi{t), then evolves further according to the full evolution 
operator U{t, r). 

Given the structure (5.11) and (5.12) of V(r), the operator J\f{ti+i,ti) has the form of 
a matrix in labels of the color basis elements, 

AA(t2,ti)|{p, /,c',c}„) 

= N{{pJU;{c,cU,W,cU;t2,ti)\{pJ,c',cU) . (^-^6) 

{c',c},n 



-20- 



Here the matrix elements are obtained from the operators h that define V(t), 



Texp 



t2 



dr h{{p,f}m,T) 



|{c}m)({c'}m|Texp 



Tdr h^{{p,f}m,r) 
Ju 



{c',c}m 



(5.17) 



Thus the Sudakov factor breaks into two factors, one for the ket state and one for the 
bra state. That is, there is Sudakov factor for each color ordered quantum amplitude. 
That is remarkably simple. However, the structure of the Sudakov factor for each quantum 
amplitude is remarkably complicated since it is an operator that mixes color basis states. 

5.8 Probability conservation 

We relate V{t) to %i{t) using the requirement that showering not change the probability 
of the hard scattering event that initiates the shower. Probability conservation requires 
eq. (2.6), 

(l|l^,(t)|{p,/,c',c}„) = (l|V(t)|{p,/,c',c}^) . (5.18) 
Here (l stands for the measurement of inclusive probability. 



c',c}m) = {{c!]m\{c]m) 



(5.19) 



The inner product ({c'}m|{c}m) corresponds to taking the trace of the color density oper- 
ator |{c}m)({c'} 

For a general statistical state \p), we measure the inclusive probability that the par- 
tons are in any configuration by using the completeness sum for the statistical states 
/, c', c}m) (as defined in ref. [4]) 



= ^^^ I d{pJ}m ^ /,c',c}m)({p,/,c',c}, 

where, to be precise, 



(5.20) 



{C,c'}r 



/ d{p, f}m=Y, 5^ / d?/b n * 

J f JO r JO : 1 

/a /b »— 

/ m 

X (27r)^<5 [pe.+Ph-'^P' 



E 



d^Pi 



(5.21) 



1=1 



Thus the inclusive probability that the partons described by \p) are in any configuration 



is 



i^P) =Y.^['^^P^f^"^ E {WU\{cU){{pJ,c,cU\p) . (5.22) 
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5.9 The inclusive splitting probability 

Using eq. (5.22), the inclusive probabiUty for a spUtting starting from \{p, f,c' ,c}m) is 

{l\ni{t)\{p,f,C,c}m)=-. \— j d{pj}^+i V {{c}„^+i\{c}m+l) 

("^ + ')'^ (5.23) 

X {{vJ ,c}m+l\Hl{t)\{p, f ,C ,c}m) ■ 

We are particularly interested in the color structure of this. We note by examining eq. (5.7) 
that on the right hand side of this equation the following generic color factors will occur: 

({c}m+l|{c}m+l)({c ,c}m+i|tjj(g)t^|{c',c}m) • 

{c,c'}m + l 

If we take the color trace of eq. (5.9) and compare it to the trace of eq. (5.10) we have 

({c'}m+l|{c}m+l)({c',c}m+i|t|4 t^|{c',c},„) = ({c'}„,|tgt^|{c},„) . ^5 ^4) 

{c,c'}m + l 

If we use eq. (5.7) in eq. (5.23) and use the identity (5.24), we find 

{l\ni{t)\{pj,c',c}m) 

= [d{pJ}rn+lS{t-T{{pJU+l)){{p,f}rn+l\ri\{p,f}rn) 



Lk 



nc{a)nc{b) riij]\^ /a/yl('?a, ^F)/b/B('7bi ^J'\) 

nc{a)nc{b) fja^T]], fa/AiVa., IJ^Dh/BiVh, ^J■F) 

1 



X - 

2 



e{k = I) 9{U+i / g) wii{{p, fU+i) (5.25) 

+ e{k = I) e{L+i = g) fU+i) - wF/'^°°"H{p, fU+i)] 

- 9{k / /) eifm+l = g)A'i^.{{p}m+l)wfl^°^''{{p, f}m+l) 

X {{c}rn\tk{fk fk+g)t]{fi fi + g) 

+ ti{fi^ fi+g)i{fk^ fk+g)\{cU) ■ 
5.10 Structure of the virtual splitting operator 

Given the structure of V{t) in eq. (5.11) and the definition (5.19) of (l|, the right hand 
side of eq. (5.18) is 

{l\V{t)\{pJ,c',cU) = HL{{pJU;{c}m,{c}m;t){{c'U\{c}m) 

(5.26) 

+ Y HR{{p,f}rn;{c'U,{c'U-,t){{c'U\{c}m) • 

Using the definition (5.12) of the color operator h, this can be written 

(l|V(t)|{p,/,c',c}™) = ({c%|/i(fe/}™,0 + /i^te/}m,t)|{c}„) . (5.27) 
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Lk 



X - 
2 



Here the operator on the right hand side is {h + i(j)) + {h"^ — icj)) = h + h'^ . Thus this relation 
determines /i + /i^ but not (f). We have assumed that </> = 0. In section 10, we will explore 
the color phase (/>. 

We demand probability conservation, eq. (5.18). Comparing eq. (5.27) to eq. (5.25), 
we see that probability conservation holds if 

K{p, /}™, = E ^ / '^^P^ f}m+iS{t - T{{p, fu+i)) {{p, fU+i \ri\{p, fU) 

nc{a)nc{b) rjs.rji^ fa/A{f]a, /^|')/;,/^(^b, A'l') 

nc(a)nc(6) r^a^ fa/AiVa., IJ'F)fb/BiVh,I^F) 

^ ^eik = I) 0(/„+i / g) wuiip, fU+i) (5.28) 

+ eik = I) 9{fm+l = 9) N/({P, /}m+l) - ^^^/'^°'^"'({P, hm+l)] 
- e{k ^ I) 0{fm+l = 9)A'i,,{{p}m+i)wfl^°^''{{p, f}m+l) 
X ti{fi fl + /m+l)4(/fc ^ fk + fm+l) ■ 

Thus we define V(i) according to eqs. (5.11) and (5.12) with h{{p, f}m,t) given by eq. (5.28). 

In the color factor, we have a factor Then is given by the same expression with 
a factor In fact, h = because tg^tj = For k = this is obvious. For k ^ we 
are dealing with a gluon with color index a exchanged between the lines k and /. We note 
that tf^tj inserts SU{(i) generator matrices in the appropriate representations on lines / 
and k, then sums over a. The generator matrices are self adjoint and they commute with 
each other because they act on different parton lines, so that = 

This result might perhaps be regarded as elegant, but it poses difficulties for practical 
implementation in a parton shower Monte Carlo program. The difficulty comes from the 
fact that the operator ti{fi — )• fi+g) t\{fk fk+g) is represented by a non-diagonal matrix 
in the standard color basis. In the full h{{p, f}m) we have a sum of such operators with 
momentum dependent coefficients. For the Sudakov factor (5.17), we need the exponential 
of an integral over t of this matrix. Thus the Sudakov factors are complicated and it is not 
easy to see what to do with them. 

5.11 Color suppression index 

To help with our analysis, we define a quantity P that we can call the color suppression 
power and a related but more useful quantity / that we can call the color suppression 
index. 

To define the color suppression power, we first note that shower evolution can generate 
a factor l/Nc" that comes from directly from applying the splitting operator to the color 
states |{c}m) and expanding in color basis states |{c}m+i)- Thus pc is the power of l/N^ 
that is associated with the quantum amplitude. 

For all of this paper except section 10, the sole source of non-zero pc is g ^ q + q 
splittings. In a 5 — >■ g + g splitting, there are two ways to connect the new q and q to the 
previous color state. The original gluon color was defined by inserting t?- along a color 3 
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line. Now with the g ^ q + q splitting we have a color matrix t°^jt1,j, in the color amplitude. 
We can use the Fierz identity, 

^ijti'j' = 2 ^ij'^i'j - ^ij^i'j' , (5.29) 

to write the result expanded in color basis states. We let pc represent the number of times 
during the shower evolution that we pick up a 1 /Nc factor by using the second term in the 
Fierz identity. 

Now we can define the color suppression power. Suppose that after enough splitting 
steps to obtain m final state partons, we reach a color density operator |{c}m)({c'}m|- 
Then we can define P{m) as the number of powers of 1/Nc in the color overlap function 
plus the sum of the explicit powers, pc, oil/Nc that are associated with the color amplitudes 
\{c}m) and ({c'} 

(^f (M„.|W™>^p{l.0(^)} , (5.30, 

where cp{m) ^ 0. The color suppression power is of obvious interest, but is less useful 
than we would like because cancellation among terms in the expansion of ({c'}m|{c}m) 
can lead to P being larger than it is for individual terms in the expansion. 

We can improve on the definition by defining a color suppression index / according to 

(i)"(M„|W„W, = |g{l + 0(-^)} . (5.31) 

where ci{m) ^ 0. The change here is that we calculate the color overlap function by using 
the color group U{Nc) instead of SU{Nc). To calculate {{c'}m\{c}m)jjrj^ y one uses the 
Fierz identity once for each gluon line, omitting the 1 /Nc term. 

The color suppression index has two properties that make it quite useful. First, as the 
number m of final state partons in the shower increases, we always have 

P{m) > I{m) . (5.32) 

That is, if we use I{m) to estimate the amount of color suppression, we can never overes- 
timate. Second, at each stage of the shower, the color suppression index either stays the 
same or else it increases: 

I{m + 1) > /(m) . (5.33) 

Thus we can think of I as measuring color disorder, like entropy: it can never decrease as 
the shower evolves. Both of eqs. (5.32) and (5.33) can be proved in a fairly straightforward 
way. We omit the proofs. 

5.12 How the color suppression index grows 

In this section, we investigate how the color suppression index I(m) grows. 

Imagine calculating {W}m\{c}m)ij(^^ y Look at any one of the gluons in the state, 
say gluon /. Its color factor was Yla^i'i^j'j- With the U{Nc) approximation, this becomes 
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(1/2) 5i'j Sj'i- Counting the rest of the graph, the indices can be connected in two possible 
ways 



In the first case, let us call gluon / a healthy gluon. In the second case, let us call gluon / 
a frail gluon. ^ With a healthy gluon, we get a factor for the sum in eq. (5.34). With a 
frail gluon, we get a factor for the sum in eq. (5.35). 
Let us look at this in more detail. 

Consider eq. (5.34) for the healthy gluon. If the healthy gluon were not there, we would 
have a factor ^ 5^' = A'c- Now insert the healthy gluon and divide by Cp oc Nc to 
normalize the states with one more gluon according to eqs. (2.4) and (2.5). The new overlap 
then has a factor N^/Nc = Nc- That is, the color suppression index has not changed. 

Now consider eq. (5.35) for the frail gluon. If the frail gluon were not there, we would 
have a factor ^ da djj = N^. Now insert the frail gluon and divide by Cp oc Nc to normalize 
the states with one more gluon. The new overlap then has a factor Nc/Nc = 1. That is, 
the color suppression index has increased by 2. 

We can now look at what happens when a frail gluon splits and what happens when a 
healthy gluon splits. 

Suppose first that gluon / is a frail gluon and splits, creating another gluon m + 1. 
Within the LC+ approximation, there are two configurations for the two gluons. The first is 
the parallel configuration, fet^"^'']*'* [^^^"li'i- With this configuration, an easy calculation 
shows that both gluons in the new state are healthy gluons. Furthermore, /(m+1) = I{m). 
The two gluons can also be in the crossed configuration, X^a \P°'^^^j'3- With this 
configuration, an easy calculation shows that both gluons in the new state are also healthy 
gluons and again I{m + 1) = I[m). Thus splitting turns frail gluons into healthy gluons 
and leaves the color suppression index unchanged. 

Suppose next that gluon / is a healthy gluon and splits. Within the LC+ approxi- 
mation, there are again two configurations for the two gluons. The first is the parallel 
configuration, ^[^"i^li'i With this configuration, an easy calculation shows that 

both gluons in the new state are healthy gluons and the color suppression index is un- 
changed: I{m + 1) = I{m). The two gluons can also be in the crossed configuration, 
Safe[^°^^]«'« With this configuration, we find that both gluons in the new state are 

now frail gluons and the color suppression index of the state increases: I(m + 1) = /(m) + 2. 

When gluon / splits, creating another gluon m+1, there will normally be other gluons 
already present and the splitting of gluon I can change the health status of some of these 
other gluons. Consider another gluon with label s. When gluon / splits to gluons / and 
m + 1 in the parallel configuration, so that /(m + 1) = /(m), a simple argument shows that 
the health status of gluon s remains the same. Now suppose that gluon / is a frail gluon 
and splits to gluons / and m + 1 in the crossed configuration, so that I[m + 1) = I{rri). 

*We call the gluon frail because if we keep the full gluon color matrices we get ^t'i'^'j'j ^ ^i'i ^j'j ~ 0- 




(5.34) 



or 




(5.35) 
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Then if s was healthy it remains healthy and if s was frail it may remain frail or it may 
become healthy, depending on how it was connected. Finally, suppose that gluon I is 
a healthy gluon and splits to gluons / and m + 1 in the crossed configuration, so that 
I{m + 1) = /(m) + 2. Then if s was frail it remains frail and if s was healthy it may remain 
healthy or it may become frail, depending on how it was connected. In summary, then, 
when a new gluon m + 1 is added so that it is color connected to gluon / in both the ket 
state and in the bra state, then splittings that increase the color suppression index may 
change healthy gluons to frail gluons while splittings that leave the color suppression index 
unchanged may change frail gluons to healthy gluons. 

We will use these observations in appendices A and B. 

6 The LC+ approximation 

With the generic structure of a (spin averaged) leading order shower set up, it is pretty 
simple to define the LC+ approximation. In the definition (5.7) of 7ii{t), in the terms that 
represent quantum interference between emitting a gluon from parton / and emitting the 
same gluon from helper parton k we replace 



in {c}m and otherwise. Thus in the LC+ approximation we keep only the color states 
in which partons / and m + 1 after the splitting are color connected. We generalize this to 
include the possibility that k = I hy replacing 



with a generalized definition of C{l,m + 1). For k = I, whenever parton m + 1 is a gluon, 
partons m + 1 and / after the splitting are always color connected, so we want C{l,m + 1) 
to return the color state unchanged. Also, for k = I whenever // = g, partons m + 1 and / 
after the splitting are also color connected, so we want C(/, m + 1) to return the color state 
unchanged.^ In the case k = I with fi = q and fm+i = Q, or vice versa, the two daughter 
partons may not be color connected. We simply define C{l,m + 1) to give one in this case. 
Thus we define 



The operator C{i,j) is a projection operator but it is not an orthogonal projection operator 
because the basis states are not orthogonal. Thus C{i,j) / C{i,j)^. 




t\{fk A + /m+l)|{c}m> ^ C(/,m + l)4(/fc ^ fk + fm+l)\{c}m) ■ 
{{c}m\tk{fk fk + fm+l) {{c}m\tk{fk ^ fk + fm+l)C\l,m+ 1) 



(6.2) 



C{i,j)\{c}m+l) 



\{c}m+i) partons i and j carry color {3, 3} or {3, 3} 
< |{c}m+i) i and j color connected in {c}m+i 
otherwise 



(6.3) 



^This case can occur for an initial state splitting with /; = fm+i ~ q with q being a quark or antiquark 



flavor. 
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Thus the LC+ approximation for T-Lj is only slightly modified from the full Tij in 
eq. (5.7): 



({p, /, c', cU+i \n'^''^{t) I {p, f, c', c}m) 

= ^Sit- T{{p, fU+i)) (m + 1) {{p, fU+i \ri\{p, fU) 



l,k 




(6.4) 



+ e{k = I) e{f^+i = g) [wu{{p, fU+i) - wt 



-eikonal 



{{Pj]m+l)] 



+ e{k / 0(/„+i = g)j^i^{{p}m+i)w\ 




{{P: /}m+l) 



) {c , c}m, ^; ^, {/}m+l) • 



Here is the color matrix defined by 



Mc({c , c}m+l, {c', C}m, U k, {/}m+l) 



= [0{k = 1)- 9{k / /)] 

X j ({c', cU+i\tl{fi ^ fi + fm+i) tkifk ^fk + L+i)CHl, m + l)|{c', c}^) (6-5) 



projection operator C'^{l,m + 1) removes the second and fourth contributions, in which 
gluon m + 1 = 3 is not color connected to parton / = 1 in the bra state. Only the first and 
third contributions remain. Similarly in figure 6 there are two contributions with full color 
but both are eliminated in the LC+ approximation. 

There are two terms in eq. (6.4) with k = I, one with fm+i / 9 and one with fm+i = 9- 
The corresponding splitting functions have collinear singularities but not soft singularities. 
In both these terms, the projection operator C{l,m + 1) acts as the unit operator, so no 
approximation is made in these terms. 

In the k / / term, fm+i = 9 and the effective splitting function A'^^wf^^"^'^ is singular 
when the gluon is soft, collinear with parton I, or both. In the limit that the gluon is 
collinear with parton I or both collinear and soft, the LC+ approximation becomes exact, 
even though these terms contain the operator C{l,m + 1). To see this, note that A'^/jufl,^"^'^ 
is independent of pk in the limit Pm+i ^Pi- Thus all of the terms with different indices 
k have the same coefficient Because of that, we can use the color identity 



+ {{c,c}m+l\C{l,m+l)tl{fk fk + fm+l)®ti{fi fl + frn+l)\{c,c}rn) 
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(6.6) 



= to write 

tkifk ^fk + g) C\l, m + l) 

k^l 

+ C(/, m + l) 4(/fc ^ /, + 5) tiifi ^h + g)] 
= -{t\Ui ^fi+g)® tiifi ^fk+g) C\l, m + l) 
+ C(/, m + l) t\{fi ^fi+g)® tiifi ^fi+g)} . 

However, after a gluon emission from line I, parton m + 1 is always color connected to 
parton m + 1 in the new color state. Thus 

tkifk -^fk + g) C\l, m + l) 

+ m + 1) tlifk ^fk+g)® tiifi -^fi+g)] ^^-^^ 
= -2tl{fi^ fi+g)^ti{fi^ fi + g) . 

That is, the operator C has no effect in the collinear or soft x collinear limit. We conclude 
that the LC+ approximation becomes exact in the collinear or soft x collinear limit. It is 
only an approximation when an emitted gluon is soft but not collinear to the emitting 
parton. 

For the virtual splitting operator, we use the the definition specified in eqs. (5.11) and 
(5.12) to write V^'^^{t) in terms of an operator h^'~"^[{p, f}m) on the space of color vectors, 
as in eq. (5.27). With full color, h{{p,f}m) was given by eq. (5.28). Now with the color 
operator C{l,m + 1), eq. (5.7) for Tif is replaced by eqs. (6.4) and (6.5). Thus we get the 
LC+ approximation for h{{p, f}m), 



'{{P, f}m, t)= yZ—, I d{p, f}m+l 
^-^ ml J 



/iLCH 

l,k 

X 6{t - T{{p, fU+l)) {{p, f}m+l\ri\{p, f}r 



nc{a)nc{b) rjarj^ /^/^(^a, MF)/b/s('?b) A^p) 



X - 
2 



nc(a)nc(6) 7?a% fa/Aiv^, ^J'F)fb/B{r|h, ^-p) 

e{k = l)e{fm+ly^ g)wil{{p,f}m+l) 

+ = I) eCfm+1 = g) [WU{{P, hm+l) - Wt°''-\{p, f}m+l)] 

- e{k / e{fm+i = g)A'i,{{pU+,)wfr'%{p, fu+i) 

X tlifl -)■ fl + /m+l) C{1, m+l) tlifk fk + /m+l) . 

(6.8) 

Now the operator h^'^'^ contains operators that act on the color space. However, the 
color basis vectors are eigenfunctions of these operators. In the case k = I, for which 
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C(/, m + 1) acts as the unit operator, the operator ti{fi — )■ fi + fm+i) i\{fi ~^ fi + fm+i) has 
eigenvalue Tr, Cp, or Ca depending on the flavors of fi and fm+i and fi. (The Ca case is 
illustrated in figure 13.) For the case k ^ I, the operator ti{fi fi+g) C{1, m + 1) tl,{fk 
fk + 5) is a Uttle more complicated. Because of the projection operator C{l,m + 1), this 
operator gives zero unless the helper parton k is color connected to the emitting parton I 
in the state |{c}m) on which h^'^'^ acts. When / and k are color connected, one gets an 
eigenvalue Ca/2 or Cp depending on the flavor of parton /. (The Ca/2 case is illustrated 
for in figure 15.) Thus 



/i^^+({p,/}™,t)|{c}™> = ALC+({p,/,c}^,t)|{c}„> 
where the eigenvalue A'^*^"'" is 
A^^+(fe/,c}™,i) 



3.9) 



y- 

^ ml 



d{p,f} 



rn+l 



l,k 

X 5{t - T{{p, fU+i)) {{p, fU+i\Vi\{p, fU) 

nc{a)nc{b) ri.^r]\^ fa/A^'Ha., l^'F)fb/B^^^^ 
nc{a)nc(b) fis,r]\, fa/AiVa., fJ'F)fb/B{'nh, fJ-p) 

X ^ 9{k = I) eifm+l + g) Wll{{p, f}m+l) 
+ 0{k = I) e{fm+l = g) [Wu{{p, f}m+l 

+ e{k / e{fra+l 
X x{k, I, {c}m) N{k, I, {f}m+l) 



(6.10) 



^-rcikonal 



g)^;,({p}„+i)wfr'^({p,/} 



({j5, f}m+l] 
m+l) 



The color eigenvalue specified in the last line is zero unless partons k and / are color 
connected in {c}m or A; = /: 



X{k,l,{c}m) 



1 k = l 

1 k and I color connected in {c}r 
otherwise 



When / and k are color connected, the eigenvalue depends on whether k 
flavors in the splitting: 



(6.11) 



/ and on the 



A^(A:,/,{/Wi) 





k = l,fi^ gjm^ 




Cf 


k = l,fi= gjm^ 




Cf 


k = lji^ gjm^ 


-1=9 


Ca 


k = l,fi= gjm^ 


-1=9 


Cf 


k^ljly^ gJmA 


-1=9 



^Ca/2 k IJi = gjm+i = g 



(6.12) 
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Since the color basis vectors are eigenvectors of /i^^"*" , the Sudakov factors for the color 
ordered amplitudes, eq. (5.17), are simply numerical factors: 



Texp 



exp 



t2 



|{c}m> 



(ir ALC+({p,/,c}„,r) 



(6.13) 



|{c}m> 



This is just the same as in the leading color approximation. In fact, with a suitable 
adjustment of what numerical factors one uses in the leading color approximation, the 
LC+ Sudakov factors are exactly the square root of the Sudakov factor for the leading 
color approximation. 

7 Weights 

Let us now look at a splitting step in the evolution equation (5.14) in some detail, using 
the LC+ approximation. With a starting state /, c', c}m) at time toj we have (using 
eq. (5.20) in the LC+ version of eq. (5.14)) 



= AALC+(t',to)|{p,/,c',c}^) 



t' 



1 



Ito (m + 1) 



d{p, f}m+l J2 ^ ^ 

l,k {c,c'}m + l 



{t',t)\{p,f,c',c}.m+l) 



(7.1) 



X ({p,/,c',cWi|7^L^+(t)AALC+(t,to)|{p,/,c',c}^) . 
Here we have explicitly displayed the sum over parton indices in T-Li: 

l,k 

The I and k dependent splitting operator is, from eq. (6.4), 

{{p, f, c', cu+i i^[:,^+(t) I {p, /, c, cu) 

= 6{t - T{{p, /}„+i)) (m + 1) {{p, fU+i \Vi\{p, fU) 

^ nc{a)nc{b) r)a.r]h /aM('7a, MF)/b/s('?b) Mf) 
nc(a)nc{b) fji^fji, /aM(^a, A^|')/f,/B(%, A^l') 

9ik = l)e{fm+i^ g)wii{{p,f}m+i) 

+ e{k = I) e{L+i = g) [mi{{p. fU+i) - ^r"({p, fU+i] 
+ 9{k + 1) eiL+i = g)A'iMpU+i)wfr"'{{p, fU+i) 

X Mc{{c,c}m+l,{c',c}m,l,k,{f}m+l,{f}m) ■ 

Here Mc is the color matrix defined in eq. (6.5). 



(7.2) 



(7.3) 
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Now consider the no-splitting operator M^'~'^{t,tQ). The color basis states are eigen- 
states of this operator, 

M'^''+{t,to)\{pJ,c',cU) 



exp 



fdr (ALC+({p, /, c}m, r) + ALC+({p, /, c'U, r)) 

■J to 



(7.4) 



X \{P,f,c',c}m) 



Using eq. (6.10), the integrand in the exponent can be written as an integral over momenta, 
a sum over flavors, and a sum over parton labels I and k: 

ALC+({p, /, c}„, r) + ALC+({p, /, c'U, r) 



^ d{pJ}ra+iY^ 



(7.5) 



X 6{T-T{{pJu+i)) {{p, fU+i\ri\{p, fU) 

X A({p, /, C}m, I, k, {p, f}m+l) ■ 

where 

A({p, /, C}m,l, k, {p, f}m+l) 

nc(a)nc(6) 7?a??b /aM(??a, IJ'F)fb/BiVh, fJ-p) 

X - 6'(A; = l)9{f,n+i / 5')uJ«/({p,/}™+i) (7.6) 
+ 0(fc = /)^(/„+i = g) [wu{{p, fU+i) - w'il^°^^\{p, /}„+i)] 

X [x{k, I, {cU) + x(A:, {c'}m)]N{k, I, {/Wi) • 

Here, x{kj,{c}m) was defined in eq. (6.11) and N{k,l,{f}m+i) was defined in eq. (6.12). 

We see that the function X{{p, f,c}mj,k,{p, f}m+i) that appears in the Sudakov 
exponent is almost the same as the integrand in the matrix element of Ti^'^'^ (t) . In fact 



{{P, f, C , C}m+l \nl^^it)\{p, f, C, CU) 

= 6{t - T{{p, /}™+i)) (m + l){{p, f}m+i\ri\{p, fU) 

X A({p, /, c}m, I, k, {p, f}m+l) C{{c, c}m+l ; {c , cjrn: ^) k, {/}„i+l) 



(7.7) 



where 



C{{c,c}m+i ){c , cjmi ^) k, {/}m,+ l) 

, {c', c\mi U k, {/}m+l) 



(7.8) 



[X{k, /, {c}m) + Xik, I, WUmik, /, {/}„+i) 



Note that Mc({c', c}m+i, W , c}m, U k, {f}m+i), eq. (6.5), in the numerator of C is nonzero 
only if one or both of x{k,l,{c}m) and xik,l,{d}m) are nonzero, so that we are never 
dividing by zero in a nonzero contribution to the sum over l,k. 
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When we insert eq. (7.7) into eq. (7.1), we see that one can generate a sphtting in 
standard Monte Carlo style by choosing the new momenta and flavors together with / and 
k with a probability proportional to X{{p, f,c}m,l,k,{p, f}m+i)- This leaves a sum over 
the choices of colors, 

^ \{P,f,C,c}m+l)C{{c,c} 
{c,c'}m+i 

For each choice of {c',c}m+i, there is a color factor C and a statistical state vector 
\{p, /, c', c}m+i) that is the input to the next splitting. In a computer program, one could 
imagine implementing the sum by summing the results returned by a splitting function 
that is called recursively. However, this is not really practical. Instead, one can perform 
the color sum Monte Carlo Style. One chooses {c', c}m+i with a probability pc, normalized 
to 

) = 1 . (7.9) 

{c',c}m + l 

Then one assigns a weight factor to the splitting equal to 

Wc{{c,c}n,+i,{c,c}rnJ,k,{f}rn+l) = — , , , T • (7.10) 

Pc[W, C|m+i, |C', C}m, h k, {f}rn+l) 

Averaged over many trials, this reproduces the desired sums. 

The shower starts with a color weight factor of Co({c', c}no) from the calculated color 
density matrix for producing an initial color configuration of Nq final state partons at the 
start of the shower. At each splitting, we multiply by the weight factor Wc from eq. (7.10). 
The shower evolves with successive splittings until it reaches some cutoff hardness value. 
At that point, let us say that we measure the expectation value of some color operator Oc- 
Then this expectation value is 

Tt[o,\{c}n){{c}n\) ={{c}n\0,\{c}n) . (7.11) 

We simply multiply the weight by this factor. In particular, making no measurement of 
color corresponds to multiplying by just the color overlap function ({c'}Ar|{c}jv). Thus if 
the final color is not measured, the complete generated event comes with a weight 

/( n Ir T \ f TT C'({c , cjm+l, {c , c}m, ^mi ^mi {/}m+l) \ / r i i n 

wtot = {{c }n\{c}n) 11 , Co({c , c}jVo) . 

(7.12) 

In eq. (7.12), the color overlap function ({c'}Ar|{c}Ar) is of crucial importance. If 
{c'}n = {c}n-, this factor is 1 or very close to 1. If {c'Iat / {cIat, this factor is proportional 
to l/ivf(^)-P^], where Pc is the exphcit power of 1/A^c that is associated with the quantum 
amplitudes |{c}Ar) and ({c'Iat] as defined in section 5.11 and P{N) is at least as big as 
the color suppression index I{N) from eq. (5.31). Recall that the color suppression index 
either stays the same or else increases at each parton splitting. 

One can use the behavior of the color overlap function to supply a hint about how 
to choose the probabilities pc- For numerical efficiency, one wants the dispersion of the 
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weights not to be large. Consequently, it is sensible to make pc small for any choice of 
color configuration {c',c}m+i that increases the color suppression index. For instance, for 
the splitting shown in figure 7, the LC+ approximation retains the first and third terms 
on the right hand side; one would pick the color configuration shown in the first term with 
probability close to 1/2 and one would pick the color configurations shown in the third 
term with probability proportional to We explore this idea further in appendix A 

using information from section 5.12 about how the color suppression index grows. 

It is also possible to choose pc = for any color choice {c', c}m+i that makes /(m + 1) 
greater than some predetermined limit /max on the color suppression index. This is an 
approximation beyond the LC+ approximation. We simply throw away configurations 
that have too many powers of We explore this idea further in appendix B 



8 The end of the shower 



At some point, the perturbative shower should end because the hardness scale encoded 
in the shower time t is not hard enough for a perturbative treatment to be reliable. At 
this point, one has a definite color density operator |{c}Ar)({c'}Ar| and one should multiply 
the color weight factor by the expectation value ({c'}Ar|Oc|{c}Ar) of whatever operator 
corresponds to the measurement to be made on the final state, according to eq. (7.11). 

In the simplest case, the operator in question is simply the unit operator and one 
multiplies by ({c'}Ar|{c}Ar). However, one may want to apply a hadronization model to 
the final state. Commonly, the hadronization model starts from the assumption that a 
color strings join the outgoing partons. A string configuration is labeled by a color label 
{cf Itv of the same form as the labels for our color basis states. Then the string state should 
be chosen with a probability proportional to the expectation value ({c'}Ar|0({cf}Ar)|{c}7v) 
of the operator 0{{ci}N) that measures whether the quantum system is in color string 
state {ci}n- 

What is this operator? The leading color guess is that it is |{cf}Af)({cf}7v|. However, 
this guess can be correct only in the leading color approximation because our color basis 
states are not exactly orthogonal to one another (and also the closed string states are not 
exactly normalized). Thus we need another set of basis states |{c}Ar,_L) that equal our 
color basis states to leading order in 1/Nc but are exactly orthonormal, 

({c}^,±|{c}7V,±) = < rn.r, ■ (8-1) 

[0 {c}n^{c}n 

The new states should be related to our basis states by a matrix equation 

|WiV>= ^^({c}7V,{c}7v)|{c}w,^) , (8.2) 

{c}jv 

where ^ is a unit matrix up to corrections. Eq. (8.1) implies that 

A\{c'}n, {c}n) Ai{c}N, {c}n) = Gi{c'}N, {c}n) , (8.3) 

{c}jv 
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where 

G{{c'}n,{c}n) = {{c}n\{c}n) ■ (8.4) 

There is a minimal solution to the equation A^A = G above. We let vl be a real symmetric 
matrix equal to 

A = VG . (8.5) 

To find A exactly, one can diagonalize G and replace each eigenvalue A by \/A. Alterna- 
tively, the expansion 

A = l + \{G-l)-l{G-lf + --- (8.6) 

allows one to write j4 as a power series in 
We conclude that it is sensible to define 

0{{ci}N) = \{cf}N,l-){{cf}N,±\ ■ (8.7) 

Then we choose the color string configuration {cf } n for hadronization with the color weight 
factor (^{c'} N\{ci} N , {{ci} N , -^\{c} n) ■ One needs to sum over the choices for color string 
configuration {c{}n. As in previous steps in the shower development, one can perform the 
sum Monte Carlo style, choosing configuration {cf}N with a probability /o({cf}iv) and 
multiplying by a weight factor 

Wc[{ci}n,{c,c}n) = 7r\\ ■ (^•^) 

P{{Ci}N) 

Note that the numerator of Wc involves matrix elements of j4. It is presumably sufficient 
to calculate Wc to leading order in l/N^- One possible string configuration is always 
{c{}n = {c}n, for which (^{cf}Ar, _L|{c}Ar) 1. Another is {cf}^ = {c'}n- These are the 
leading possibilities if ({c'}7v|{c}Af) oc 1/A^c with p = or 2. For p = 4, there can be other 
important configurations such that both ({cf}Ar, _L|{c}Ar) oc l/N^ and ({c'}Ar|{cf}Ar, _L) oc 
1/Nl 

One does not necessarily need to follow the LC+ shower all the way to a (1 GeV)^ 
virtuality scale. One could use the LC+ shower for a few splitting steps, down to, say, a 
(100 GeV)^ scale. Then one could assign a classical color string configuration {cf}]\f to the 
partonic state as outlined above. After that, one could use a leading color shower to get 
from the (100 GeV)^ virtuality scale to the (1 GeV)^ scale. This would be appropriate if 
the measurement to be made on the final state is only minimally sensitive to what happens 
at the softer scales. 

9 Full color perturbatively 

The shower evolution equation with full color has the form 

\p{t)) =l({t,to)\p{to)) , (9.1) 
where U{t,to) obeys the evolution equation eq. (2.2), 

j^U{t,to) = [ni{t)-V{t)]U{t,to) . (9.2) 
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We have approximated U{t,to) by U^'^^{t,to) where 

jU^^+{t,to) = - VL^+(t)]Z^LC+(t,to) . (9.3) 

This differential equation can be solved iteratively in the form eq. (5.14) 



Here AA'^*-'+(t2, ti) is the no-splitting operator, 



, ti) = exp 



(9.5) 



It is well to recall here the essential point: the operator V^'~^~^{t) is diagonal in the standard 
color basis that we use, so that it is practical to calculate its exponential. 
Now, what if we want shower evolution with full color? Then we need 

j^m.to) = [?^LC+(i) _ vLC+(^) + ^^^(^) _ AV(0]^(t,to) , (9.6) 

where 

/LC 



(9.7) 



(9.8) 



AV(t) = V(t)-VL^+(t) . 
This evolution equation is equivalent to 

U{t,to)=U'^'^+{t,to)+ fdTU{t,T) [Ani{T)-AV{T)]U'^'^+{T,to) , 

Jto 

which can be solved iteratively: 

Jto 

+ [ 'dT2 Hdn U^'^+iti, T2) {AUi(t2) - AV(t2)] ^/^^+(r2, Ti) 

X {AUiir^) - AV{t^)\U^'^^(t^M) 
+ ••• . 

One can have any (small) maximum number of insertions of [A'H7(r) — AV(r)] . For 
instance, to start with one might test whether one such insertion makes a significant dif- 
ference. To make one insertion of \AUi{t{) - AV(ti)] , one would generate t\ at random. 
Then one would run an LC+ shower from the starting scale to up to scale t\. After that, 
application of AHiijx) produces a weight and a new shower state, which is the starting 
point for an LC+ shower from t\ to the final shower time tf. Similarly, application of 
AV(ri) produces a weight and a new shower state, which is the starting point for a second 



(9.9) 
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LC+ shower from ri to if. The resulting values for the measurement function from the 
two second stage showers would then be summed. 

It is interesting to note the counting of logarithms in eq. (9.9). Suppose that the 
parton shower is used to calculate an observable O in which there is a large logarithm L, so 

that (O) has the form c{n, 2n)a"L^" + c{n, 2n — l)a'^L'^'^~^ H . Suppose further 

that a shower with full color, generated by U{tf,tQ), correctly calculates all coefficients 
c(n, 2n) and c(n, 2n — 1). Then the first term in eq. (9.9) correctly generates the coefficients 
c{n,2n) since the LC+ approximation is exact with respect to color for the soft x collinear 
singularities. In the second term in eq. (9.9), the one insertion of [A?^/(r2) - AV(r2)] 
generates a factor asL by correcting the color content of a wide angle soft splitting. This 
factor multiplies factors a^~^L^"~^. Thus the second term makes contributions to the 
coefficients c(n, 2n — 1). The third and higher terms in eq. (9.9) contribute to coefficients 
c{n,j) with j < 2n — 2 only. Thus to get the coefficients c(n, 2n — 1), one needs only the 
ffist two terms in eq. (9.9). 

In general, by using one or two terms beyond the LC+ approximation in eq. (9.9), one 
can see whether the splitting operators beyond the LC+ approximation have an impor- 
tant influence on whatever observable is being investigated using the parton shower. One 
may expect that for many observables, the operators AT-Li{t) or AV(t) are not important. 
Including some factors of [AT-Li{t) — AV(r)] can test this hypothesis. 

10 Soft gluon exchange phase 

Up until now, we have presented the evolution equations for a parton shower in which the 
virtual splitting operator V is determined by the real splitting operator T-lf together with 
an assumption about the structure of V. 

However, the structure assigned to V{t) leaves out an important physical effect: the 
wave function of two partons emerging from a hard interaction can accumulate an SU(3) 
phase factor ?7 ~ 1 + \(f){t)dt by exchanging a soft gluon. Here (j) is an operator on the 
quantum color state. In an inclusive cross section, this phase cancels between the virtual 
graphs for the bra state and for the ket state. That is, the phase cancels in a completely 
inclusive measurement because 



Tr 



U\{c]m){{c]m\uA = {{c]m\U^U\{c]m) = {{c ] m\l\{c] ^m) ■ (10.1) 



However, the color matrix (j){t) changes the color state, which influences the further evolu- 
tion of the shower, so that the effects of the color phase do not cancel from all observables 
measured at the end of the shower. These effects can be important [11-15]. 

Recall the structure of eqs. (5.11) and (5.12) for V{t), which we can summarize in a 
shorthand notation as 

V{t) ^ [h{{p, f}m. t) + \ct>{{p, /}„, t)]®l + l® [/it({p, /}^, t) - i0({p, t)] . (10.2) 



^"The phase is sometimes called the Coulomb gluon phase by analogy with the Coulomb phase in non- 
relativistic quantum scattering. However, in a relativistic calculation in Coulomb gauge, only part of 
the phase comes from the Coulomb force between colored partons; the rest comes from the exchange of 
physically polarized gluons. 
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The probability conservation equation determines h through^ ^ 

(l|'W/(i)|te/,c',c}„) = ({c'}^|M{P,/}m,t)+/i^({p,/},n,t)|{c},„> ■ (10.3) 

This gave us h from 'Hj{t), as in section 5.10. That is, we could find h without looking 
at any details of virtual graphs by simply knowing that the singularities of real emission 
graphs have to cancel with corresponding singularities of virtual exchange graphs. 

Previously, we made the assumption that (p = 0. However, when a soft gluon is 
exchanged between partons / and k, the graph has a part proportional to i times a hermitian 
matrix (p in the color space. This color phase is not zero. A simple calculation using the 
eikonal approximation, similar to other calculations in the literature (for example, ref. [11]), 
gives 



cPi{pJU,t)= -2n^e{k^l)"-^^^^f^n-T, . (10.4) 

rrf 47r ^ ' 

Lk 



In order to make the connection to parton showers clear, we present our calculation of this 
result in appendix C. In eq. (10.4), we sum over partons / and over partons k with k ^ I. 
Thus each pair of partons appears twice in the sum. Not all combinations contribute: 



1 k and / are final state partons 

1 k and I are initial state partons . (10.5) 
otherwise 



The phase depends on the relative velocities of partons k and I: 



Vki = \n -( ^2 — • 10.6 

V [Pk-pir 

Note that if either parton is massless, Vki = 1. There is a factor evaluated at a scale 
/XR that depends on the shower time and potentially depends on the momenta of partons 
k and /, depending the physical meaning of the shower time used in the parton shower 
(see appendix C). Finally, there is a color operator • T;. The operator inserts a color 
matrix on line k: if \I' is the vector representing the color state {c}m before the virtual 
exchange as in eq. (2.4) or eq. (2.5), then the effect of the is to map 

^j,ai--afc--a„ _^ j,b ^ ^ai---a^--a„ ^ (10.7) 

Here b is the color index of the exchanged gluon. Similarly, Ti inserts a color matrix 
on line I. Finally, we sum over the color index b, as indicated by the dot in ■ Ti. The 
operator ■ Ti maps |{c}m) into a linear combination of color basis states |{c}m). 

Let us denote the part of V{t) that comes from (j){t) by V0(t). Imagine calculating the 
expectation value of some observable O by using a parton shower in which we calculate 



^^Recall that h = with full color but that in the LC+ approximation we defined /i^*"^^ with [/i^*^^]^ 7^ 
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perturbatively in powers of V0(i). We note immediately that contributions proportional 
to an odd number of factors of V0(t) do not contribute to (O) because (O) is real and 
these contributions have an odd number of powers of i. However with an even number of 
powers of V(/,(t) we have a factor ib(i7r)^"' = ib(— l)"7r^"' and have a generally nonvanishing 
contribution. How big these contributions are depends on how sensitive the observable is 
to the color flow of the event. This is a question that deserves further study that is beyond 
the scope of the present paper. However, one can note immediately that the development 
of a parton shower depends crucially on the color structure of the shower state because 
this color structure determines the preferred emission direction for relatively soft gluons. 
Thus a sudden change in color structure caused by inserting two occurrences of V(f,{t) at 
some early shower time t can have a substantial influence on the flow of momentum at the 
end of the shower. For instance, a gap in the rapidity of emitted partons could be created 
or could be filled in. 

Can we say more about the color structure of the phase operator (/)? We find that the 
color operator T^. • T; applied to |{c}m) gives 

Tk ■ Ti\{c}m) = k, {cU) xik, I, {cU)\{cU) 

(10.8) 

+ ^ 0, k, {CU)\{CU) + ^ T^\{c}m)j . 

The first term does not change the color state and does not change the color suppression 
index. In the first term, xik, I, {c}m) is one if / and k are color connected, zero otherwise, 
as in eq. (6.11). The eigenvalue A is 

{2 k,l carry colors {8,8} connected to a two 
parton closed string . (10.9) 

—1 otherwise 

The second term is a color suppressed contribution, with a factor l/N^ and a factor 
C{l,k,{c}m) = 



1 k, I carry colors {3, 3} or {3, 3} 
-1 /c,/ carry colors {3,3} or {3,3} . (10.10) 
otherwise 



The corresponding contribution to the gluon exchange phase cancels exactly between the 
phase 4>{{p, f}m,t)\{c}m) of the amplitude and the phase —{{c'}m\4'{{p, f}m,t) of the 
conjugate amplitude. The remaining term, T^, changes the state |{c}m) to a linear com- 
bination of other color basis states. 

The contribution to from the first term in • T/ in eq. (10.8) can be included in 
the LC+ approximation V'^*^'^ for V. Then the remaining part of becomes part of AV 
and can be treated perturbatively as in eq. (9.8). 



'^^We count the 1/A'^c factor that multiphes as increasing the color suppression power pc by 1 in 
eqs. (5.30) and (5.31). With this definition, the term {l/Nc)T^i either increases the color suppression 
index, eq. (5.31), of the color state {c,c'}m or leaves it unchanged. 
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The color phase operator (f) depends on the parton masses. There are two places 
where the mass dependence might be important for the analysis of TeV scale processes at 
the LHC. First, in the later stages of showering, one can produce gluons with transverse 
momenta not too far above the 6-quark mass. Sometimes such a gluon can split to 6 + 6. 
Then the h and h are non-relativistic, so the mass dependence matters. Second, the hard 
process under investigation can produce top quarks or perhaps squarks, gluinos, and other 
very massive particles. In these cases, the particle masses matter. 

Even though particle masses can matter, it is of interest to understand what happens 
when all parton masses can be neglected, so that v/fe — > 1. There can still be dependence 
on the parton labels A;, Z if the relation between the shower time and the renormalization 
scale in a.^ depends on the parton kinematics. Let us suppose that we neglect any such 
dependence. Then 



/}m, *) = - 2vr g ^ 7^ /) Oki Tk ■ Ti 



(10.11) 



Since the whole shower is invariant under color rotations, we have used color basis states 
that are overall color singlets. Applied to color singlet states, we have 







(10.12) 



for any I. Prom this, we derive 



(f>{{P^ f}m, t) = 0o({P, f}m, t) + (l)a.hi{p, f}m, t) 



(10.13) 



where 



^o({p, f}m,t) = 2tt 



47r 



as 



kh{{pj}m,t)= -2tt ^4Ta-Tb 



(10.14) 



In the first term in 00, we sum over final state partons while in the next two terms we sum 
over the initial state partons. In either case, we have a color operator Ti-Ti, which is just 
the Casimir operator, with eigenvalue Ca or Cp depending on whether parton i is a quark 
or a gluon. Thus 



)({P, f]m,t)\{c}m) 



as 



Nl) + C^{N^^-Nl)\\{c}m) , (10.15) 



where A^J is the number of gluons in the final state of {p, /}„j while is the number 
of gluons in the initial state and is the number of quarks and antiquarks in the final 
state while is the number of quarks and antiquarks in the initial state. When we apply 
—(pQ to the bra state ({c'}to|, we get exactly the opposite phase. Thus the term c^o in 
the color phase contributes nothing to the development of the shower and can be simply 
dropped. This leaves the single contribution ^abj representing an effective double strength 
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color exchange between the initial state partons. Notice that with the approximation used 
here there is no color phase for e^e~ annihilation or for deeply inelastic scattering. 

Using eq. (10.8) for the color operator Ta-Tb, we obtain a leading color term that leaves 
the color state unchanged plus a term that changes the color state. The leading color term 
provides a phase that occurs if partons a and b are color connected in the state {c}m- There 
is an exactly opposite phase that occurs if partons a and b are color connected in the state 
{c'}m- Thus there is no net phase if partons a and b are color connected both in state {c}m 
and in state {c'}m and there is no net phase if partons a and b are not color connected 
in state {c}m or in state {c'}m- If partons a and b are color connected in one of {c}m or 
{c'}m but not the other, then there is a net phase that appears in every evolution interval 
until the color connection situation changes. These phase factors will tend to reduce the 
contribution of this sort of state to whatever observable is to be measured. 

There are also contributions to Ta • that change the color configuration. These terms 
have the potential to change the energy fiow in the final state by changing the color flow. 

11 Conclusions 

In general, a parton shower Monte Carlo event generator should generate contributions to 
the density operator in color space in which the color in the ket vector |{c}m) and the 
color in the bra vector ({c'}m| are different. Standard parton shower generators based on 
evolution from hard splittings to soft splittings work in the leading color (LC) approxima- 
tion, which is the leading order in an expansion in powers of ^/N^- In the leading color 
approximation, only states with {c'}m = {c}m occur. 

In this paper, we have introduced the LC-|- approximation, a generalization of the 
leading color approximation. Going beyond the leading color approximation inevitably 
involves sums over color choices. These sums can performed with Monte Carlo summation: 
selecting a choice at random according to a prescribed probability function. At the end of 
the shower each event comes with a weight. 

The LC+ approximation has several nice features: 

• For each splitting, the leading soft x collinear singularity and the leading collinear 
singularity are treated exactly with respect to color. There is an approximation with 
respect to color, but it occurs only in wide angle soft splittings. 

• Evolution can start with any state |{c', c}m) . Thus if one starts with a hard scattering 
process, one can fully use the color subamplitudes that multiply color states |{c}m) 
for the hard scattering, including in the calculation interference between ket states 
|{c}m) and bra states ({c'}m| with difi^erent color configurations. 

• The Sudakov factors are numbers. That is to say, the standard color basis states 
|{c', c}m) are eigenstates of the Sudakov operator. One does not have to exponentiate 
non-diagonal matrices in the color space. 

• In fact, not only are the Sudakov factors numbers, but also there is a separate Sudakov 
factor for the ket state |{c}m) and for the bra state ({c'}m|. This feature may prove 
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useful for matrix element matching when the color structure of the amplitudes is 
treated exactly. 

• With a simple extension of the formalism, one can include in the LC+ approximation 
the phase induced by exchange of soft gluons. 

• The LC+ approximation is still approximate: remainder terms in the generators of 
shower evolution are left over. However, the remainder terms can be included in a 
perturbative calculation up to some order. 

• The LC+ approximation can provide an efficient tool to sum large logarithms with 
full color at leading and next-to-leading log level for a certain class of observables if 
one uses the first perturbative correction as described in section 9. 

The inclusion of weights in the shower generation has the potential to produce numeri- 
cal problems. If the dispersion in the weights is large, then the number of events needed to 
calculate the expectation value of some observable with reasonable accuracy will be large. 
The total weight for an event is the product an initial weight, a final weight, and weights 
for individual splittings. By multiplying a large number of individual weights, one has the 
potential to produce large weights. For instance, (l + l/N^)'^ can be large if is large. For 
this reason, it seems likely that there is a practical maximum to the number of splittings 
that can be generated using the LC+ approximation. 

Fortunately, it is possible to turn the LC+ approximation off before it goes too far. 
In fact, we have found two ways to do that. First, one can simply run the LC+ algorithm 
for some number A'max of splittings and then return to the LC approximation. For that, 
one must replace mixed states |{c',c}m) by diagonal states |{c, c}m). We have presented 
a plausible model for doing that. Second, one can continue using the LC+ approximation, 
but not allow the generation of states with more than a certain amount of color suppression 
as measured by what we have called the color suppression index. 

We are thus encouraged that the LC+ approximation can prove useful. The authors 
do not have at immediate hand a dipole based Monte Carlo event generator suitable for 
implementing the approximation.^^ However, the LC+ approximation is quite general and 
could, we think, be implemented in an existing generator. 

Finally, one can ask if one really needs to go beyond the LC approximation. The 
answer is that we do not know if one needs to go beyond the LC approximation. The LC 
approximation is the first term in an expansion in powers of We do not know if 

the contributions of order or beyond are important. One could imagine that these 

contributions are not important for some observables but are important for others. If 
that is the case, we do not know which observables are sensitive to l/N^ effects. An 
implementation of the LC+ approximation would give us the means to investigate these 
questions. 

^^We are, however, working on such a program. 
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A Choice of probabilities 

As outlined in section 7, a practical shower evolution computer algorithm implementing 
the LC+ approximation can use the same method as with a leading color shower to chose 
the index I of the parton that splits and the index k of the helper parton as well as the 
momenta and flavors of the daughter partons in the splitting. It remains to chose the color 
configuration {c',c}m+i of the daughter partons. Here, standard methods do not suffice 
because the functions that play the role of quantum probabilities in the leading color 
case are now not always positive. One can, however, still use the Monte Carlo method 
for implementing the sums over colors, choosing the new color configurations at random. 
The program should chose color configurations {c',c}m+i at each splitting according to a 
probability function Pc{{c' , c}m+i, W , c}m, l, k, {f}m+i)- This function is a probability, so 
it must be non-negative and normalized to 

X] Pc{{c,c}m+l,{c,c}m,l,k,{f}m+l) = l- (A.l) 

{c',c}m + l 

Then, following eq. (7.10), at each splitting we multiply the accumulated weight for the 
shower history by a weight factor Wc for the splitting given by Wc = C/pc, where C is the 
color factor given in eq. (7.8). 

There is an initial weight factor Co({c', cjiyo) from the calculated color density operator 
for producing an initial color configuration for A'^g final state partons at the start of the 
shower. 

At the end of the shower with final state partons, there is an final weight factor 
corresponding to the observable to be measured. If the observable is not sensitive to color, 
then the part of this final weight factor that is associated with color is the color overlap 
function ({c'}Ar|{c}Ar). 

Thus the complete generated event comes with a weight as given in eq. (7.12), 

/( n Ir T \ f TT C'({c , cjm+l, {c , c}m) ^m, ^m; {/}m+l) \ / r i i \ 

wtot = {{c}n\{c}n) 11 f Co({c,c}jVo) , 

\m=I Pc\l'^ 1 C/m+li iC , C\mi I'm-, >^rm \J jm+lj / 

where the color factor C is given in eq. (7.8). The expectation value of the observable 
calculated after many shower events is independent of the probability function pc used 
in the calculation. However, the dispersion of results from the many shower events does 
depend on the choice of pc- In order to have a low dispersion of values and thus an efficient 
calculation, one wants the net weights Wtot not to vary too much. In particular, the weights 
Wtot should never be too large. 
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In order to keep Wtot from being large, it seems that the factors pc in the denominator 
of eq. (7.12) should be large, but there is a constraint from the normalization condition 
eq. (A.l): the probability budget should not be spent on configurations where it is not 
really needed. For that reason, one would set Pc = for color configurations {c', c}m+i for 
which C = 0. 

One can also make pc small for color configurations that will, in the end, make 
{W}n\{c}n) smah. Now, {{c'}n\{c}n) oc (l/A^c)^('")-^^ where Pc is the amplitude 
power defined in section 5.11 and P{m) is the color suppression power defined by this 
relation. Recall that P{m) > I{m), where I{m) is the color suppression index defined in 
eq. (5.31). The color suppression index either stays the same or grows at each splitting 
step, so it is sensible to make pc small for any choice of color configuration {c', c}m+i that 
makes / grow. 

We consider a g ^ g + g splitting for the case that k ^ I and that parton k is also a 
gluon. Look at Mc, the numerator of C({c', c}m+i, {c' , c}m, Im, km, {f}m+i), as defined in 
eq. (6.5): 

Mc{{c,c}m+1 , {c , c}m) I, k, {/}m+l) 

= -({c , c},n+i 14(5 ^g + g)^ tk{g ^g + g)C\l, m + i)|{c', c}^) (A.2) 
- ({c',c}™+i|C(/,m + 1)4(5 ^ g + 9)(^ti{g ^ 9 + g)\{c!,c}m) . 

Recall that C(/,m + l) selects color states in which partons / and m + 1 are color connected. 
The operators t\ and t\ here insert the new gluon into the color states according to 

t\{g^g + g) = ^/c;{al{l)-al{l)) . (A.3) 

Here a!^(/) inserts gluon m + 1 just to the right of gluon I on its color string in |{c}m) and 
aL(Z) inserts gluon m + 1 just to the left.^^ This notation is explained in more detail in 
sections 7.2 and 7.3 of ref. [4]. The factor y/Up comes from the normalization of the color 
states, eqs. (2.4) and (2.5). Thus 

Mc{{c,c}m+i , {c , cjmj I, k, {/}m+l) 

= -Cf{{c' ,cU+i\{al{l) - alii)) {a+{k) - a-{k))CHl,m + l)\{c' ,cU) (A.4) 
- CF{{c',c}m+i\C{l,m + l) (alik) - al(A:)) ® (a+(/) - a-{l))\{c' ,c}m) • 

It will be useful to define 

if k lies just to the right of / on a string in {c}m 
otherwise 

(A.5) 

if k lies just to the left of / on a string in {c}m 
otherwise 

^^Here "left" and "right" refer to the listing of gluon attachments in eqs. (2.4) and (2.5). The adjoint 
operators a±{l) act analogously on the states ({c'}m,|. In drawing a picture to represent the bra color state 
({c'}m|, one normally reverses the direction of the arrow on the color 3 line. However, this does not change 
the definition of "left" and "right." 



X+ikJ,{c}m) 
X-{k,l,{c}m) 
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Then 

{{c}m\a±{k)C\l,m + l) = {{c}.,n\a^{l)XT{k,l,{c}m) , ^^^^ 
C{l,m + l)aiik)\{c}m) =XTif^,l,{c}m)a^:^il)\{c}.m) ■ 
Using this notation, we have 

; {c , c}m; ^) k, {/}m+l) 
= ({C ,c}rri+l|a^(/) a+(/)|{c,c}m)CF[x+(^,/,{c}m) + X+ ^ {c }m)] 

+ ({c',cWi|al(/) ®a_(0|{c',c}™)C7F[x-(A;,/,{c}„) +X-(fc,^{c'}m)] (A.7) 

- {{c,c}m+i\a^+{l) ^ a^{l)\{c ,c}m)CF[x+{k,l,{c}m) +X-{k,l,W}m)] 

- ({c ,c}m+i|al(0 ^ a+{l)\{c ,c}m)CF[x-{k,l,{c}m) +x+{k,l,{c}m)] ■ 

What does eq. (A.7) teh us? There are functions x±{k, I, {c}m) and x±{k-, I, {c}m) that 
describe the color structure of the initial state. One of the functions x±ik, I, {c}m) equals 1 
if parton k is connected to parton / in the state |{c}m). One of the functions x±{k, I, {c'}m) 
equals 1 if parton k is connected to parton I in the state ({c'}m|- Potentially, parton k 
is color connected to parton I in both the bra and the ket state. If parton k is not color 
connected to parton / in either the bra and the ket state, then Mc = 0. Then there is a 
color factor C-p. Finally, there is a factor that describes the color structure of the state 
after the splitting. For instance, ({c', c}m+i |a^(/) a+{l)\{c' ,c}m) is 1 if gluon m + 1 is 
inserted to the right of gluon I in the ket state and in the bra state and this factor is zero 
otherwise. There are four possible ways to insert the new gluon, corresponding to the four 
terms in eq. (A.7). The first two terms correspond to parallel configurations while the last 
two terms correspond to crossed configurations. The fact that the crossed configurations 
come with minus signs tells us that we cannot use the coefficients in eq. (A.7) as Monte 
Carlo probabilities for the respective insertions. 

We recommend defining the probabilities pc according to the status of gluon I in the 
color configuration {c', c}m- If gluon / is a healthy gluon, then the color suppression index 
remains unchanged if we chose one of the parallel configurations for the insertion of the 
new gluon, but if we chose one of the crossed configurations, then the color suppression 
index increases by 2. Thus we can chose parallel configurations with a probability close to 
1 times an overall normalization factor and we can choose crossed configurations with a 
probability proportional to times the overall normalization factor. We define 

Pc{{c',c}m+i ) {c , c}m, I, k, {/}m+l) 

= ({c ,c}m+i|a^(/) (g)a+(/)|{c',c}r„)— — [x+(A;,/,{c}m) +x+{k,l,{c}m)] 



+ 



2Cf 



{{c', C}m+l\ai{l) (g) a_(/)|{c', C}m) ^^-^ [X-{k, I, {c}m) + X-{k, I, {c'}m)] (-^ g) 



CaA 



+ ({c',c},„+i|a;^(/) ® a_(0|{c',c}m)^^-^[x+(fc,Z,{c}m) + X-(/i:,^,{c'}m)] 
+ {{c ,c}m+i\aL{l) a+{l)\{c' ,c}m) ^ '^^2 [x-ik,l,{c}ni) + X+{k,l,{c'}m)]; 
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where the normahzing factor \p is 

\p = [x{k, /, {c}^) + x{K I, {c'}m)] ■ (A.9) 

Here x{k,l,{c}m), defined in eq. (6.11), is 

X{k,l,{c}m) = X+ik,l,{c}m) + X-ik,l,{c}m) ■ (A.IO) 

This is the counting factor that appears in the denominator of C in eq. (7.8). 

If gluon / is a frail gluon, then the color suppression index remains unchanged for 
both parallel and crossed insertions of the new gluons. Therefore, we allow all of the color 
choices with equal probabilities: 

Pci{c,c}m+i , {c , c}m,, k, {/}m+l) 

= ({c,c}m+i|a;^(/) (g)a+(/)|{c',c}m)-'^[x+(/!:,/,{c}m) + X+(^,Z,{c'}m)] 

+ {{c\c}rn+l\al{l) (g) a-{l)\{c\c}rn)y[X-{kJ,{c}rn) + X-{k,l,{c'}rn)] -|^-|^^ 

+ ({c ,c}m+i|aV(0 ® a_(0|{c',c}m) — [x+(fc,/,{c}m) + X- (^, ^, {c }m)] 

+ ({c,c}m+i|al(0 ® a+(0|{c',c}m) — [x-(fc,/,{c}m) + x+(fc,^{c}m)] . 

Here the normalizing factor is 

X'p = 2[x{k,l,{c}m) + xik,l,{c'U)] . (A.12) 

B Avoiding too high color suppression 

The LC+ approximation has the good feature that it can generate a shower starting from 
the color density operator terms |{c}m)({c'}m| with {c'}m 7^ {c}m- Such a state inevitably 
leads to {c'Iat 7^ {c}n at the end of the shower and thus to a weight factor ({c'}iv|{c}Ar) 
proportional to to some power. However, the LC+ approximation can generate color 

states with overlaps proportional to to a very large power. Suppose that we are 

content to limit the accuracy of shower evolution (within the LC+ approximation) to 1 /Nc 
to a certain power /max- For instance, one might set /max = 4. Then we can expand 
weight factors ({c'}m|{c}-m) as a power series in 1/A^c and not calculate contributions with 
too many powers of 1/Nc- Furthermore, one can track the color suppression index /(m) of 
generated states. One can arrange not to generate any splitting that makes /(m+1) > /max- 
Let us see how to avoid generating states with I{m + 1) > /max- Suppose that we 
have a color state with /(m) = /max- There may well be some frail gluons in the state. 
Within the LC+ approximation, these gluons can split. As we have seen in section 5.12, 
the splitting of the frail gluons leaves /(m) = /max and turns the frail gluons into healthy 
gluons. 

There will likely also be some healthy gluons in the state. As we have seen, when 
a healthy gluon splits in the parallel configuration, the color suppression index is left at 
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I{m + 1) 

— -^max- However, when a healthy gluon sphts in the crossed configuration, the 
color suppression index increases to /(m + 1) = /max + 2. This is beyond the limit that we 
want to maintain. Thus, when I{m) = /max; the splitting of a healthy gluon in the crossed 
configuration should not be allowed. 

This analysis suggests that one should modify the probability function pc for a healthy 
gluon splitting when I{m) = /max by dropping the last two terms in eq. (A. 8), so that one 
now uses 

, {c , c}m) h k, {/}m+l) 

= ({c , c}m+i\a^+{l) (E) a+(/)|{c', c}m) Y [x+{k, I, {c}m) + X+{k, I, {c}m)] 
+ {{c',c}m+i\aL{l) (g) a-{l)\{c ,c}m)^[x-{kj,{c}m) + X-{k,l,W}m)] ■ 

The corresponding normalizing factor Ap is 

~Xp = xikJ,{cU) + xik,l,{c}m) ■ (B.2) 

This leads to a weig ht factor = 2Cf/Ca = l-l/N^. Since this procedure in any case 
drops color suppressed contributions once /(m) has reached /max, it makes sense to simply 
change the weight factor to = 1. 

If one were to set the probability pc to generate a splitting of a healthy gluon in the 
crossed configuration to pc = € <^ 1, then one would very rarely generate an event with 
I{N) > /max with a very large weight. We set e = so that no such events are actually 
generated. One should note, however, that then there is no control over contributions 
to the calculated observable at order greater than l/AT^™""". Thus, for instance, one can 
calculate the color overlap ({c'}Ar|{c}Ar) at the end of the shower as a power series in 1/Nc 
and drop terms of order greater than l/A'^c'"'"'"^''. (See eq. (5.31) for pc-) 



C Calculation of gluon exchange phase 

To define the phase from soft gluon exchange, we consider a gluon with momentum q = 
{E, q) exchanged between possibly massive final state partons / and k. Using the eikonal 
approximation, and supplying an ultraviolet cutoff mi and an infrared cutoff m2, the 
amplitude for this process is 

(2vr)4 J J [pi-q + le) [pk ■ q - le) (^ + le) 

Here pi and pk are the momenta of partons / and k as they enter the final state. We let 
pi = Eivi and pk = EkVk where 

Vk = (l,Wfc) . 
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Thus vi and Vk are the 3-velocities of the particles in the reference frame that we are using. 
If either of particles has a non-zero mass, then the absolute value of the corresponding 
velocity is less than 1. Thus our integral is 

(27r)3 y V {E-vi-q + ie){E-Vk-q-ie){E^-q^+ie) ^ ' 



Now we can perform the E integration to give 

" (2^)2 J <m< ^i) 



{vk-vi)-q + ie\\^+vi-q \^ - Vk ■ < 

We are interested in the imaginary part of G, which comes from replacing 1 / [{vk —vi)-q + ie\ 
by -ITT 5{{vk - vi) ■ q): 

nu,G= -i^^^2TrTk e{m2 < \q\ < mi) 

X 6{{vk -vi) -q)! ^ — =, + 



\^+Vi-q \^-Vk-q 

Now choose coordinates with both Vk and vi in the x-z plane, with z in the direction 
of Vk — vi and vi^x > 0, Vk^x > 0. Then with q± = (q^, Qy) and Qx = —\q± \ cos 6*, we have 

4TT\vk-vi\ I^-lI Jo l-vi^xCosO 

In eq. (C.5), there are two terms that are related by g — )• —q and vi o Vk, but since 
vi,x = Vk,x, these two terms are equal. We have just written one of the terms here and 
multiplied by 2. 

The integrals are trivial: 

ilmG = - 27ri "l^' 2T; • Tfc log(mi/m2) . (c.7) 

47r\vk - vi\[l - vfj^/^ ^ ' 

The velocity dependent prefactor can be simplified using 

\vk-vi'^\"^-v\^ = {vk-vifvli , (C.8) 

where 

' ,,2 ,,2 / ^2 ^2 



V (^^fc-^'O V iPk-pi)^ 



Then 



iImG= -2^i-^2Tr Tfc log(mi/m2) . (C.IO) 
4vr Vki 

If the soft gluon is exchanged between initial state partons / and k, there are some sign 
differences in the derivation but the the result is the same. However, if the soft gluon is 
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exchanged between an initial state parton and a final state parton then we find ImG = 0. 
Thus the general result is 

iImG= -2m— ^2r;-Tfc logfmi/ma) , (C.ll) 
Vkl 47r 

where 6ki is 1 if A: and I are both final state partons or both initial state partons and zero 
otherwise, as in eq. (10.5). 

The logarithm in eq. (C.ll) results from an integration between m2 and mi. 

If we say that \q\ is related to shower time t by 

1^1 =CQoe-* , (C.12) 

where Qo represents the hard momentum scale corresponding to shower time zero and C 
is any dimensionless constant, then = —dt, so 

iImG= - 27ri — ^ 2T; • Tfc dt , (C.13) 
Vkl 47r Jt^ 

where mi = C Qoe~^^ and m2 = C Qo^"*^ • The coefficient of dt here is the differential 
gluon exchange phase, 

i^iki{pJU,t)= -2m^^2Tk-Ti , (C.14) 

Vkl 47r 

as reported in eq. (10.4).^^ 

We have related the shower time t to |^ in eq. (C.12) with the aid of an unspecified 
parameter C. Different parton shower schemes use different choices for shower time. For 
our present purposes, we do not need to make a definite choice. Nevertheless, it is useful 
to present at least one choice in order to illustrate the physics issues involved. 

Note that there is an energy denominator for the parton I plus gluon state before 
parton / absorbs the gluon, previously emitted from parton k: to first order in \q\/Ei, we 
have 

AE = ^pf +mf- \q\ - ^{pi + q)^ + mf ~ -{\q\ + Vi ■ q) . (C.15) 

This is the denominator in the first term in eq. (C.5). One sensible way to define the 
shower time for the first term in eq. (C.5) would be to set 

\q\+vi-q^QQe-' . (C.16) 

That is, integrating over a range of shower time up to t2 would mean integrating over q 
subject to\(^+vi- q> Qo exp(— ^2)- Using the change of variables in eq. (C.6), the relation 
between shower time and the integration variables is 

\q\{l-vi^^cose)^Qoe-' . (C.17) 

^^Note that in eq. (10.4) the term corresponding to gluon exchange between partons I and k occurs 
twice in the sum over parton indices. That is why an explicit factor 2 disappears between eq. (C.14) and 
eq. (10.4). 
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It is convenient to define 9 by 

vi^^ = cose . (C.18) 
(For massless partons, 6 is half the angle between vi and v^-) Then the relation is 

2\q\\sm^{e/2)+cose sm'^{e/2)\^Qoe-^ . (C.19) 

We integrate over 9. If 9 is of order 1, then the integration over 9 is dominated by 9 values 
of order 1. If <C 1, then the integration over 9 is dominated by 9 values of order 9. Thus, 
after the 9 integration, the effective relation between \q\ and the shower time is 

2|glsin2(^/2) ~ Qo e~* . (C.20) 

Thus one possible relation between t and \q\ is eq. (C.12), with C = l/[2 sin^(0/2)]. 

What about the renormalization scale, //^(t)? One might base this on the denomina- 
tors of covariant Feynman diagrams instead of the energy denominators. Define Di(q) = 
(Pl + q)'^ ~ ^fi where pi is the final state momentum of parton /, pi = (Ei,ps) with 

— [Pi + rnf\''^l'^. The exchanged gluon is virtual, but let us consider the on-shell 
contribution with E = Then 

Di{q) = 2pi-q = -2Ei{\c[\+Vi-q) . (C.21) 

If we use —Di for we would take //^ = 2£^;Qo exp(— t). Similarly, we might take 
/x^ = 2SfcQo exp(— i). For the graph as a whole, 

^il = {El + Ek)Qoe-' (C.22) 

might seem sensible as long as Ei and E^ are not too different. 
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